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SUMMARY 


A study is presented of the technical feasibility of growing single- 
crystal silicon ribbon in space environment. Electromagnetic shaping and 
other shaping processes are considered in separate tasks Procedures are 
described for calculating the electromagnetic fields produced in a silicon 
ribbon by an rf shaping coil. From a knowledge of these fields, the forces 
on the ribbon and the degree of shaping to be expected are determined. fne 
expected 9teady-state temperature distribution in the ribbon is calculated 
in the one-dimensional approximation. It is indicated that shaping coil 
current required has in the past frequently been underestimated, currents 
in 100 A range being needed to produce ribbon of 0.4 mm thickness, even with 
very small shaping coils very close to the melt. The power requirements for 
shaping an isolated ribbon are rather large, by the standards of near-term 
space missions, but for shaping to 0.4 mm a ribbon being drawn from a free- 
floating spherical melt, the power required is in the range of a few kilo- 
watts which does not seem excessive. Calculations on simplified models 
indicate, however, that lack of flatness of the shaped ribbon and excessive 
heating of the melt by the eddy currents induced by the shaping fields may 
pose problems An analysis of the relative effects of various kinds of forces 
other than electromagnetic showed that in the space environment capillarity 
forces would dominate, and that the shape of the melt is thus principally 
determined by the shape of any solids with which it comes in contact This 
suggests that ribbon may be produced simply by drawing between parallel vires 
A concept has been developed for a process of off-angle growth, in which the 
ribbon is pulled at an angle to the solidification front. If it could be 
implemented, such a process promises to offer increased growth rate, better 
homogeneity, and thinner ribbon. While we recommend that work on development 
of concepts for space processing of silicon ribbon be continued to the extent 
consistent with the overall goals of the space-processing program, we do not 
recommend further terrestrial experimental work because of the highly space- 
oriented nature of the most promising concepts. 
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1. INTRODUCTION 


Background 


From almost every aspect, the least satisfactory part of the process 
of manufacturing silicon single-crystal wafers, whether for integrated-circuit 
or solar-cell application, is the cutting and polishing of the wafers from 
the crystal boules. The process is expensive, requires skilled hand work, 
wastes at least half the crystal, and introduces strain and heavy-roetal 
impurities into the finished wafers. Consequently, there has been consider- 
able interest in developing methods of producing single-crystal ribbons or 
films directly. In this report, we describe investigations into some aspects 
of whether processing in a space environment might be a feasible and useful 
method of producing silicon single-crystal ribbon. A little interest also 
attaches to other semiconductors, but inasmuch as silicon is by far the most 
widely used material for Integrated circuits and is the prime candidate for 
photovoltaic cells for large-scale power conversion, it will be the only 
material specifically discussed. Some other work has been done along these 
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lines. A program at Texas Instruments^ looked into electromagnetic shaping 

of silicon ribbon; it involved some terrestrial experiments which were not 

( 2 ) 

too successful. A more recent program at McDonnell-Douglas is devoted to 
a conceptual and economic study of an overall space processing system. In 
this proposed system, electromagnetic shaping is again used, but the 
ribbon is drawn from a continually repelenished levitated sphere. This con- 
cept is very nicely adapted to space processing. We shall describe aspects 
of these programs in some more detail later. 

1.2 Organization of Program and of Report 

This is the final report on a program involving three separate 
tasks, the first two of which were carried out separately. The first task 
was a study of electromagnetic shaping, including development of methods of 
calculation of the electromagnetic fields, forces, and power requirements, 
and thermal effects involved. This work is described in Section 2 , where 
we show among other things that the current needed for electromagnetic shaping 
has generally been rather underestimated in previous work. In the second 
task, we examined other possible space-processing methods and isolated the 
principal factors involved in space growth of single-crystal ribbon. The 
results of this investigation and a suggested method for drawing ribbon in 
space are described in Section 3. The third task, discussed in Section 4, 
was an assessment of the overall technical prospects for space production of 
silicon ribbon and a recommendation of the directions future work might 
profitably take. 

Several areas were specifically excluded from consideration because 
it was felt investigation of them would be premature: 

1. mechanical stability against dynamic perturbations, 

2 . degree of crystallinity and perfection of ribbon produced ^except for 
a few general comments), 

3. economic feasibility of any process. 
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In this report, equations and figures are numbered separately 
for each section, using a multiple decimal system, in order to avoid 
excessive renumbering. Thus, Eq. (2.3.1) refers to the first equation 
of section 2.3. The references, on the other hand, are numbered consecutively 
for the whole report to avoid excessive repetition. Since there are only a 
few tables, they too are numbered consecutively through the whole report. 

1.3 Physical Properties of Silicon 

In Table I we have collected, fox convenience, representative data 
on the physical properties of silicon which are potentially useful in the 
study of crystallization and electromagnetic shaping in this material. Not 
all the data in Table 1 are essential to the calculations presented in this 
report, though. We have not attempted to make a complete literature survey 
or to determine a "best" value for the parameters, but we have insofar as 
possible selected experimental results determined near the melting temperature 
in both the solid and liquid phases. These data will be entirely accurate 
enough for the present studies. Here, and throughout the report, SI units 
will be used. 

We shall adopt the following convention for the symbols referring 
to these data; symbols without a subscript will refer to the liquid phase; 
symbols referring to the solid phase will have a subscript "s". 

The thermal dif fusivities , x * K/CA, calculated from me Table I 

-"5 2 —6 2 

data are h * 3.0 x 10 m /sec In the liquid state and K g * 9.9 x 10™ m /sec 

for the solid. The latter value can be compared with an extrapolated estimate 

—6 2 ( 13*i 

of 9.4 x 10 in /sec given by Touloukian et al^ . These dif fusivities are 

large, considering the small dimensions over which shaping of a ribbon is 

liable to take place; so in investigating thermal properties, it will not be 

necessary to consider other :han steady-state conditions. 

The surface tension, and consequently power required for shaping, 

might be reduced by heating the melt above the melting temperature, particularly 

( 2 ) 

reduced pressure . We will not consider this possibility in our numerical 
calculations, though. 



4 


TABLE I . TYPICAL VALUES OF PHYSICAL PROPERTIES OF 
SOLID AND LIQUID SILICON NEAR THE MELTING TEMPERATURE 


Proparty and Symbol 

Value 

Refs. 


Liquid 

Solid 


Melting temperature, 

1685 1 

K 

3 

Latent heat ot fusion, L 

1.80 MJ kg” 1 

3 

Density, A 

2.49 x 10 3 kg nf 3 

2.29 x 10 3 kg m" 3 

4 

Thermal conductivity, K 

67. W rtf 1 K _1 

22.0 W m" 1 K _1 

5,6 

Specific heat, C 

9.12 x 10 2 J kg -1 K _1 

9.73 x 10 J kg" K' 1 

3 

Resistivity, p 

8.3 x 10 ^ ohm m 

1.7 x 10 3 ohm m 

7 

Surface tension, y 
(or surface free energy) 

0.72 N rtf 1 

1.23 N m" *(f 111} plane'i 

8,9 

Viscosity, |i 

2.0 x 10 3 N sec m ^ 


10 


Emissivity, e 


0.5 (rough estimate) 0.6 (extrapolation) 


11,12 
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2, ELECTROMAGNETIC SHAPING 
2.1 Statement of Problem 

The objectives of this portion of the program are the development 
of analytical techniques and ca lculational procedures, including computer 
programs, for determining electromagnetic fields and forces and temperature 
changes within a column of molten silicon being shaped into a thin ribbon by 
eddy-current forces induced in the silicon by an external shaping coil; 
determining the power requirements and degree of shaping that may be obtained 
for given configurations and dimensions of the shaping coil: and assessing 
effects of the freezing oi the shaped ribbon near the shaping position. 
Electromagnetic shaping is obviously attractive as a potentially contamination- 
free shaping method, and appears particularly attractive for use in space, 
where there will be no difficulties with convection in the melt. Moreover, 
it has generally been considered, apparently on the basis of the early 
estimates of Gaule' and Pastore^^, that the power requirements for electro- 
magnetic shaping to submilliroeter thickness were quite moderate. If the 

ribbon were drawn in space from a large levitated sphere, as in th: ^^onnell- 

( 2 ) 

Douglas concept , further savings in required power might be anticipated 
because of the removal of the constraint of constant cross-sectional area of 
the molten zone which the gravitational field effectively imposes on earth. 

One might in principle attempt to solve the electromagnetic problem 
in a completely general way for a given coil configuration by assuming 
some reasonable shape for the molten ribbon, calculating the eddy currents 
and force distribution for this shape, allowing the shape of the melt to be 
modified slightly in the direction indicated by the force distribution, 
recalculating the currents and forces, and so forth, until the increase in 
surface free energy of the melt on further shaping Just balances the electro- 
magnetic work necessary to achieve that shaping. Actually we ahall follow a 
program rather like this, but to avoid having to solve an extremely complicated 
three-dimensional eddy-current problem, in conjunction with a free variational 
problem to determine the shape ot the melt, we will make what we believe are 
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reasonabTe assumptions about the geometry of both the shaper and the ribbon. 

This will also enable us to identify mure readily the factors of primary 
importance in determining the effectiveness of the shaping process. In 
return, we give up in principle, though not in practice, since the exact 
problem would surely prove wholly intractable, some knowledge about the 
exact shape of the frozen ribbon. 

As far as the shaping coil goes, it is clear that it is sufficient 
to consider it as made up of pairs of straight wires running parallel to the 
flattened surfaces of the ribbon and perpendicular to the pulling direction, 
and symmetrically disposed across the melt, since wires running in any other 
directions will not be effective in producing ribbon of the desired shape. 
Moreover, it will be adequate to consider one pair of such wires, since the 
field distribution from a group of shaping wires can be obtained by adding 
vect iially the fields from each constituent. Thus we shall largely 
concentrate our attention on a single hairpin-shaped coil, with the bend 
sufficiently remote from the ribbon to have negligible influence, disposed 
as described above (and illustrated in Figs. 2.3.1, 2.3.2, and 2.4.1). 

We will neglect any effects of the induced fields back on the shaper, which 
after all carries a very high current. The single hairpin shaping 'oil is 
also very important as a limiting case; for a given total input current and 
a given distance b of closest approach of coil to melt, it is easy to see 
that no greater degree of shaping can be obtained than can be gotten by putting 
all the current into single coll at the distance b. Some more complex 
configuration may have some other desirable feature, such as shaping effectively 
over a greater ribbon length, but for a given ribbon thickness, it will 
take more current. 

As for the ribbon, we shall assume that it is drawn in *'nch a way 
that in the critical shaping region its brond sides are reasonably flat and 
’’vertical" — that is, perpendicular to the plane of the shaping coil. Calcu- 
lations on surface energy of melt-ribbon conf igurat ions v ** ' ^ seem to 
indicate that such a configuration can be obtained with a shaping coil similar 
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to that described above. Of course the coil must then be interacting 

significantly with the non-"' ro rticar’ portion of the liquid, but we must ignore 

this in the interests of simplicity. It does seen/^ that there are circumstances 

where this is not at ail a bad approximation. Another approximation is that 

any possible change of ribbon shape between the plane of the shaper, where 

the electromagnetic force is greatest, and the freezing plane is ignored. 

It might be possible under some circumstances to arrange the pull rate so 

that the flat ribbon shape was . aintained, but one would not wish to *ely 

on this. Finally, we assume the ribbon is very broad in comparison to its 

( 2 ) 

thickness, as it appears it will have to be for satisfactory productivity; 
so it will be a good approximation to assume that the fields, forces, etc., 
do not vary along the broad dimension of the ribbon. This is a tremendous 
simplification, since it reduces a fully three-dimensional problem to a two- 
dimensional one, (There may still be field components in all directions, but 
they do not vary along one). To our knowledge, no three-dimensional eddy 
current problem of comparable complexity has ever been solved, numerically 
or otherwise. The flat ribbon would presumably be terminated by little 
caps of semicircular cross-section (Fig. 2.4.1); although these may be 
important in determining the surface energy, the shaping field effects on 
them can be neglected. Originally we attempted to avoid the flat ribbon 
assumption by assuming the ribbon cross section was that cf an elongated 
ellipse. Eventually, however, we realized that this only complicated the 
problem without giving much in return; in particular, it is muon harder to 
consider the effects of the freezing interface if the ribbon is not flat. 

The effects of the various approximations are discussed qualitatively, and 
to some extent quantitatively, elsewhere in the report. 

When specific ribbon dimensions are needed for purposes of numerical 
calculation, we shall assume a ribbon width of 6 cm and thickness 400 am. The 

( 2 ) 

same values are generally used for illustration in the McDonnell-Douglas report . 
For solar cell applications, somewhat thinner ribbon-say 250 |im — would be 
de^irrthle; in the ribbon-from- levitated-sphere setup, though, the amount of 
e t i a power required to achieve the lesser thickness does not ordinarily appear 


to '>0 very great. 
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2.2 Basic Electromagnetic Equations 

In this section we give, without derivations, the basic equations 
for determining the eddy current density in and induced force distribution 
on a metallic melt in the field of external shaping coils. SI (mksA) units 
are used. A single exciting frequency will be assumed; generalization is 
straightforward. A superscript 0 will be used to indicate values of various 
qu .atitities in the n vacuum M region outside the molten zone; no superscript 
will indicate values inside the melt. We discuss first the molten ribbon 
and the space around it. When we come to discuss the solidified ribbon and 
the solid - liquid interface, we will use a superscript "s" for fields and 
forces within the solid region. A tilde (~) underneath a letter will be 
used to designate a vector quantity. 

For a single driving frequency f = w/2 rr, and neglecting, as usual 
in eddy-current problems, small displacement current terms, the independent 
Maxwell equations to be solved are: 


Inside 

Outside 


V x El =-iu uuH 
~ 0 ~ 

0 0 
v X E = - i4 Q wH 

(2. 

V x H = a E 

0 

V x H =0 

(2. 


V • J2° = 0 

(2. 

where a is the conductivity of the molten material, and jE and H denote the 
steady parts of the electric and magnetic fields; that is, E(t)= Ee icUt ,for 
example. The vsual procedure of taking the curl of the curl equations leads 
to the vector Helmholtz equations (inside) and vector Laplace equations 
(outside") 

^ 2 E = (2i/6 2 )E 
'“ w 

V 2 E° = 0 

(2. 

? 2 H = (2i/6 2 )H 

7 2 h° = 0 

(2. 

1 ^2 

where the skin depth 6 = (2/^.qHc) ' . 




2. la,b) 
2 . 2a ,b) 
2.3) 


2 .4a ,b) 
2 . 5a, b) 


where the 
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In the solid, the corresponding Maxwell and Helmholtz equations will 
naturally be of the same form as (2.2. la, 2a, 4a, 5a) , but with o replaced by 
iT g , the solid conductivity, and 6 by 4 g , the skin depth in the solid. 

The skii* depths 4 and 4 g are among the most important physical 
parameters of our problem (or of any eddy-current problem) . To provide an 
idea of orders of magnitude, we note the numerical relationship 

4.737f = 10 6 6 2 , (2.2.6) 

where f is in MHz and 4 in j,m. Thus for example for an operating frequency 
of 21.1 MHz, the skin depth in the melt will be 100 p,m. At the same frequency, 
the skin depth 4 g in the solid will be about 455 gjn. The boundary conditions 
to be satisfied at the melt-vacuum interface are 


For the electric field 

E = E ° 

•St ~t 

(V x E) - (V x E°) 

~ t ~ t 

E « E ° 
mi ~n 


For the magnetic field 

a - h° 

~-t M: 

(V X H°) = 0 


H - H 

*M1 'Ml 


(2.2, 7a, b) 

(2.2.8a,b) 
(2.2. 9a f b) 


where subscripts t and n denote as usual components tangential and normal 
to the surface, respectively. Not all these conditions are independent. At 
the solid- liquid boundary, the same equations hold if superscript f, 0 M is 
replaced by "s" . Additional boundary conditions, at large distances from the 
ribbon and at the ribbon’s center, are discussed where they are specifically 
needed . 

The eddy current density in the melt may be determined from Ohm’s 

law 

j, * a E , (2.2.10) 

and the steady force density on the melt found from 


f 



(j x H*) 


( 2 . 2 . 11 ) 
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If we define a complex Poynting vector by 

£ = 1/2 (E x H*) , (2.2.12) 

and tiien define 

n - - S • n dA , (2.2.13) 

where the integral is taken over the surface of the ribbon and n denotes 
an outward normal to that surface, then Re II is the steady Ohmic power loss 
in the ribbon and Im II is 2 uj times the average magnetic energy stored in the 
ribbon^"* J . At the input of the shaper coil, the effect of the eddy currents 
is thus to provide an additional impedance * 211 /I , where I is the peak 
current in the coil. 

2.3 Fields and Forces in Uniform Ribbon 

In this section we investigate the forces exerted on a semi 
infinite region of a uniform metallic material (such as liquid Si) by a single 
long straight wire carrying alternating current of peak value I. The wire is 
parallel to the surface of the metallic half-space, which we shall refer to as 
the melt, although a real molten metal would of course be deformed away from 
the wire. Thus, the principal assumption in addition to those made previously 
is that the liquid-solid interface is far enough away from the shaping region 
that it can be ignored. Techniques to avoid having to make this assumption 
are described in Section 2.6, and its region of validity is discussed in 
Section 2.7. Also, since we are really primarily interested in a thin ribbon 
ratnei. than a semi-infinite region, the results of this section will only be 
strictly valid when the skin depth is small compared to the plate thickness 
t say 5 <T t/3. This approximation will also he removed in Section 9 .6. As 
we Midi l see, there may be conditions under which a relatively large skin 
depth M.nv be desirable; but the results of this section will still be good 
first opproximations at any reasonable operating frequencies. 

Aspects of this problem have been studied before, particularly by 
l' 'i i ( bkv and Jerrard^*^ , and by Stoll^*^. The force distributions were not 
studied, thougn. Our simple method for obtaining the fields and forces at 
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the surface of the melt is also new. Furthermore, Stoll's account contains 
a good d«V of incorrect mathematics. We have tried to clarify this part of 
the d scmu ion. This relatively simple problem is very instructive in show- 
ing the g neral nature of the force distributions and the relevant combinations 
of tarcmer&rs of Interest. 

The notation to be used is Illustrated in Fig. 2.3.1. We call the 
wlre-tn-melt perpendicular distance b. The direction along which the melt 
ie ^eing d awn is x; this is made more apparent in Fig. 2.3.2, where the two- 
wire thin ibbon case is illustrated. By the assumed symmetry of the problem, 
nc quantii ites can vary with z. The problem to be solved is therefore two- 
dimtnslona l . 

We also note that the magnetic field of the Isolated straight wire 
has no component along the direction of current flow, (the z direction) and 
that the presence of the molten zone, with no variation along z either, can- 
not alter this. Thus our Maxwell equation 7 x H - o E inside the melt reduces 
to fc[B v /dx - BlWBy] = aE. So E inside the melt and, by continuity, outside 
has onj_y a z-component . We can write the Maxwell equations in component 
form a: 


- 1 ? 

(2.3.1a) 

BE- 

i u 0 -H y 

(2.3.1b) 

cifljj _ dH* J = aE z inside 

(2.3.2a) 

d:-c <*y i = 0 outside 

(2.3.2b) 


Rather than u?ing these equations and the associated boundary conditions 
directly, it will be cc venient to follow Stoll and introduce the magnetic 
vector potential A, iefined by B - m. q H = 7 x A. For H not to vary in the z 
direction, the x- and y- components of A must be constants, which can be taken 
without lor' of generality to be zero. Thus, we are left with 


11 

x 


1! 

V 


= 1 £Az 

By 

1 ^ 

" "u c lx 


-• - i'uA 

z 


and Lrom (2.3.1a) 


(2.3. 3a) 
(2.3.3b) 


F 


(2.3.3c) 
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y 



FIG. 2.3.1. PERSPECTIVE OF LONG STRAIGHT WIRE ABOVE SEMI-INFINITE 
MELT, SHOWING COORDINATE SYSTEM USED, 
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where again an unimportant constant term has been dropped. Differentiating (3a) 
with respect to y, (3b) with respect to x, and using (2) and (3c) we get the 
Helmholtz equations for A^, 

* q < 2 A 3 inside the melt, (2.3.4a) 


a 2 A 3 + d 2 A 3 
3x 3y2 


oA. d*A 

a 2 + *2 = 0 outside, (2.3.4b) 

where we have indicated by A^ the value of A inside the melt (region (3) in 

Fig. 2.3.2) and by A ? the corresponding value just outside. In (2.3.4a), 

2 2 1 

a = 21/6". The boundary conditions at the surface of the melt (y = 0) can 

— ♦ •'■4 

readily be obtained from those for E and H. Continuity of yields 

A 2 = A 3 , y • 0 . (2 .3.5a' 1 

Continuity of gives 

^2 = ^3 , y * 0 (2.3.5b) 

3y By 


Continuity of provides no additional information. A z (x,y) must equal 

A (~x , y) in order that H (x,y) = H (-x,y) and H (x,y) = - H y (-x,y) , as in 
z x x y j 

the case of the isolated wire. Thus also 3A /3x = 0 at x = 0. 

z 

We will solve Eqs, (2.3.4), together with boundary conditions (2.3.5) 

md appropriate conditions at infinity, by the Fourier transform method. 

>*efore doing this, however, it will be helpful to discuss the vector potential 

of the isolated current-carrying straight wire. 

_( i s ) 

The vector potential A at (x,y) due to an isolated wire at (0,b) 
is given by 


A (i8) (x,y) 


In (x 2 + 

4 r 


(y-b) ) K 


where K is an arbitrary constant vector, independent of x and y. It is 
easily shown that this yields the usual formulas for the components of the 



r**>i " ii "• 


magnetic field of the isolated wire. As with the electrostatic potential, 
a constant term may be added since only differences in potential have 
physical significance. The difference in vector potential between two 
points of different y but the same x is clearly 

A (ls) (x, yi )- A (l8) (x,y 2 )= (2. 

1 4 tt xiCyj-b)* 

/ 13 ) 

Consultation of a table of integrals (for instance, Eq. 1.4(2) of Erdelyi v ) 
shows that (2.3.6) can also be written in the form 

* <l *>(x., x ) - A (1,) (x.y 2 ) .!£?>£! co.kx<. k<r 1 - k> -, k<y 2- b) ) 
for y\ and y 2 < b. This integral is perfectly well behaved, but the separate 
exp mential terms both diverge at the lower limit. However, we will not 
get into difficulty using such integrals to represent the vector potential 
of the isolated wire at a given point--that is, 


718/ v ti 0 Ik > dk i 

* <x ’ y) ta • 


k(y-b) 


so long 8S we remember that only potential differences are physically 
significant . 

Since A is an even function of x, a cosine transform on this 
variable is appropriate. We define 


A i (k,y) 


^ A^ ( x , y) coskx dx, i - 2,3 


Multiplying Eq. (2.3.4) by cos kx and integrating on x, we find 

*. . 


r A3 iX l 2 X ~ 

£ = (k 1 Of )A. 


(2.3.7b) 
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Thus In general, 


», . C 2 « k » . 0 2 « 


and 


V C 3 e 


V k 2 + n _Vk 2 + Qf Z y , 

7 + e 


( 2 . 


( 2 . 


where it must be borne in mind that the coefficients C^, C^, D 2> can 
depend on k. In order that the solution be finite as y — ®, we must 

have * 0. The Fourier Inversion Theorem thus gives 

A 2 (x,y)*^ J Q [C 2 (k)e ky + D 2 (k)e“ ky ] cos kx dk , (2. 

2 «> >1 lr ^ 4. Qf ^ v 

A 3 (x,y) » ^ , q Cj (k) e y cos kx dk (2. 

When c r -• 0-that is, when the skin depth 6 becomes infinite — these 
results must go over to those for the isolated wire; that is 

U I *oo lf V 

A i "* 2 n~ J 0 k cos kx dk. «-* 0, i- 2,3. 

This can be achieved If 

D 2 - 0, a - 0 

C, - C., a - 0; and 
u I e- 1 * 

V k) -rV • 

C 2 (k) must be independent of a for the vector potential to approach that for 
the isolated wire at large b. Inserting (2.3.9) into (2.3.5), we see 


C 2 + °2 " C 3 


k(C 2 -D 2 ) - Jk 2 + or 2 C 3 


3.8a) 

3.8b) 

3.9a) 

3.9b) 
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Thus, ve finally obtain 

(2.3.10a) 
(2.3.10b) 


A 2 < x ,y) - IT Jo L' 


‘‘o 1 r. -, k( y- b) e‘ k(y+b) 


k — OT-a bj co * te dk * 


Hi 


k + a y-kb 


A 3 (x,y) - IT" J 0 l ~ J k 2 


dk. 


From Fqs . (2.3.10), all the fields forces, etc., of Interest can be calculated 
if the integrals can be evaluated. The integrands of A 2 and A^ are finite 
and integrable everywhere in the range of integration. 

From (2.3.3) and (2.3.10b) we can find the field components 
Inside the metallic region 

-i^Ipa, jp + y-kb 


9 n Jc 


cos kx dk 


<k‘ + AT + k 


U « - I 

X TT j 0 


i •» gT?Z + c».b.dk . 


.k 2 + Of 2 + k 


(2.3.11) 


(1.3. 12) 


H -if" e^~ kb k lin kx dk 
y tt J 0 \" W 1 "it 

,k + cr + k 




By multiplying the integrand by 'fr "i 1 7 , the alternative forma 

ik* + « - k 


(2.3.13) 


«.-£K - * <*- 177?, 


„ . «! J* * * 2 T-» fa (k 2 . .2 . »J?T7)dk, 


X 2tt1 j 0 


H + ^ Y_kb sin kx (k 2 - klk 2 + # 2 ) dk (2.3. 13') 

y 2ni .JO « 


a 


(2.3.11 ) 
(2.3.12') 


can be obtained. 



i 

I 

I 
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These forms are useful when b/6 is large. Numerical calculations^* 6,1 ^ 
show that the fields fall off very rapidly (faster than exp <-|y|/6)) Inside 
the melt. The force distribution, depending on a product of electric and 
magnetic fields, falls off more rapidly still. Thus it seems appropriate to 
concentrate our attention on the fields and forces at the surface of the melt, 
where they are largest; however, by (2.2.13), these fields are all that are 
needed to calculate the power requirements. 


At the surface, we have 


-kb 


h <°> 



cos kx dk 


I f® e ^ coskx \!k^ + a „ 

", (0) ■ 5 Jo IT— 1 * 




w + a + k 


coskx dk - 


T™ e ^ cos kx k dk | 

Jo jTT^rf : — J 

Ik + cr + k 


-V w ■ i Jo frryj 


lnkx dk 
k 


(2.3.14) 


(2.3.15) 

(2.3.16) 


where E (0) is short for E (x,y»0) , etc. These integrals can be evaluated 
z z 

in terms of certain generalizations of Bessel functions known as Struve 
(19^ 

functions , Bessel functions of the second kind, and elementary functions 
by writing the cos kx and sin kx terms in exponential form and using the 
known results^ 16 ^ 


r. 


dk e 


-pk 


k + \lk + a 


5 F“ 2 ^ % (ap) ' Y i (ap)) 


1 

2 2 

i p 


(2.3.17) 


and 


'*> -bk 

;o e 


+ X 


> 


coskx dk* 


(2.3.18) 
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l 


together with the result, obtained from (2.3,17) by differentiating with 
respect to the parameter p. 


dk k e 

Jo 


*pk 




* a ^ 


\i»p <v ap) - V ap)) • 3 -“rif • < 2 - 3 

• p 


Eqs. (2.3.17) and (2.3.19) are valid tor Re p ■> 0, |arg a| < n/2, conditions 
which are always fulfilled in the problem of interest here. In these 
equations, 11. and H. are the Struve functions of the first and second orders. 
Tbcsy should not be confused with Henkel functions which are also frequenty 
denoted by H's but without the tilde. (Also, although we used the tilde else- 
where to denote vectors, the Struve functions are of course not vectors). 

We find that the fields at the surface can be expressed a 3 


and 


11 <ul . « »2 * . i 

h <°> • -W cji 7^2 - r lx <&i< x *> - V x >» * w*-'-*,'’ 

V°> - U ii* - { f [ fe h VV <■ 

" r * } \ 

1 hi 3 s 2 -i n ,, , 

3 8 T 7b?> 3 )J • ( 

»^ }] . «-3. 


where 

- eid-?) + i(i+5)] . 

X_ - B[(i+§) ^ i(l-§>1 
e - b/6 

and 

5 * x/b 


(2.3 

(2.3 

(2.3 

(2.3 


.19) 


,»]} 

20 ) 

21 ) 

22 ) 

. 23s) 
.23b) 

.23c) 

.23d) 
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The variation In the complex variables \+, as g Increases from zero Is 
sketched In Fig, 2.3.3. It should be notad that we have expressed each field 
In terms of the product of a coefficient which Is sellar to the field for 
the Isolated wire and a function ot the dimensionless variables 0 and g. 

Thus the most appropriate combinations of variables for studying the problem 
have been Identified. Also, we may make use of the known properties of the 
Struve and Bessel functions to obtain In a straightforward way useful 
numerical results. 

We could now write down general expressions for the force distribution 
at the surface, but we find it more convenient to consider first the fields In 
the large and small 0 regions. 

For large 0 (0 » 1), X+- and X- will always be large in magnitude, 

(19) 

and we may use the asymptotic expressions 

x r r (tn 4 1/2) 

V z) - V z) • n r (n -m ♦ 1/2) (z/2) 2m ~ n + 1 * 

These expansions are valid on the entire first sheet of the complex plane 
except for the negative real axis, which, as Fig, 2.3.J shows, Is never crossed. 
Alternatively, we can set y * 0 in Eqs. (2.3.11 / -13 / ) , and expand the remaining 
square roots. Either way, we obtain the asymptotic expansions 


E ' <0> ' n oM‘-(rT{7y[ l -'i up- * w (TT^)2 ’ 1 7^)2'] , <2 - 3 - 24) 

H ,<°> -4t^L 1 -'T8 iq* ♦ 1 Qs (feSfv] • < 2 - 3 - 25 ’ 


and 


H y (0) L 1 - i 0- | T7??)] 


(2.3.26) 


The next term in H y (0) is of order 0 


-3 



-u* imrrWHH&ti&P 
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FIG. 2.3.3. THE COMPLEX \ -PLANE CUT ALONG THE NEGATIVE REAL 
AXIS, SHOWING THE VARIATION OF THE FUNCTIONS 
\ r , X_, (EQS. 2. 3. 23a, b) WITH THE PARAMETER f; . 

FOR f LESS THAN ZERO, \ + and \_ ARE INTERCHANGED. 
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For small ft— that Is, for the vire very close to the melt — we can 
use expansions' 19,20 ^ of the H' s and Y's valid for small values of the 
argument to obtain 


*,<»> 

H,(0) 


o 

2tt 

X. | 

nb - 


1 + 3 e * 1 C ln 8 J 1 ?" + v. Y - i/J 

2(1+? 2 ) + 3 f 1 + 


+ 


H (0) 

y 


i 

2TTb 



h 


In this case we have only displayed terms through 0; the first neglected 

2 

terms are proportional to B . Through terms in ft, H^(0) is real. 

Now by using the relation (Eq. (2.2.11)) 


- - - 
f - ^ Re ( j x H*1 

where 

j = a E 

we may find the force distribution at the surface in the two situations 

(large 8 and small 8) for which we have given expressions for the field 
components. Since 


(2 

(2 

(2 


E = E k and K* = H* £ + H* ^ 
z x y 

the force distribution has in general two components, one normal to the 
surface 


£ „ ■ ¥ " «,<> 


and one tangential to the surface 


U- a 

f t - 4- Re 


(2 


3.27) 

3.28) 

3.29) 


.3.30) 


2 


(2.3.31) 
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(2.3.24-26), 


(2 
(2 

There are several general observations that can be made about these force 
distributions: 

1. f^, as would be expected, is a force on the melt, tending to deform It 
away from the wire. 

2. i t is a stretching force, tending to pull the surface of the melt away 
from the line of closest approach. 

3. Both force components fall off rapidly when |^| becomes large. 

-1 -2 

4. The coefficients of 0 , B » etc., are bounded for large 5; so it is 

clear the expansion is valid for large B for all values of F . 

5. The maximum value of |f fc | occurs approximately at |f | ■ 1//7; 

the maximum of f is of course at § * 0. 
n 

6. For large g, the ratio of the tangential force to the normal force is 
very small. To be specific. 


In the asymptotic (g » 1) regime, we find from Eqs. 
(2.3.30) and (2.3.31), that 


-v 2 

1 

1 . X (13?. 

TT 2 b 2 6 

(1+? 2 ) 2 


s i 

[ i . 1 1-3C 2 ' 

ttV 

(1+5 2 ) 4 ' 

1 2 * !+S 2 ■ 


f t (0) ■“ 343/7 

.. — _ 

f (0) max 4096 

n 


ss 0.2216 

2 — Q^* 


so f max is less than 1 percent of f max for g * 5. 

t n 2 2 

In Fig. (2.3.4), we have plotted ^ n (°)/(UQ I Mb O a * a function 

of the normalized distance along the pulling direction £ for g * 5 and 

g = Ic is obvious from the curves that this quantity depends only very 

weakly on B for g £ 5, and tabulation or plotting for other B values in the 

2 2 4 

asymptotic range is unnecessary. In Fig. (2.3.5), we plot f t (0)/(^ 0 I /tt b ) 


3.32) 

3.33) 



* !(H r/ TrW) 




FIG. 2.3. j. TANGENTIAL FORCE DENSITY AT MELT SURFACE, FOR TWO VALUES 
OF 0 = b/4 , AS A FUNCTION OF F * x/b, THE NORMALIZED 
DISTANCE ALONG THE hELT AWAY FROM THE LINE OF CLOSEST 
APPROACH (F = 0) OF THE SHAPING WIRE. CURVES ARE ANTI- 
SYMMETRICAL ABOUT ^ = 0. 
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vs. F, again for 6-5 and 0 = <=. The same weak dependence on 0 is observed. 
While we have normalized the curves so as to plot quantities that are finite 
for large B, actually f - » and f -* 0 as 0 - » for fixed b. In Fig. (2.3.6), 
we plat on the same axes f n (0) /(u Q I 2 /n^b 3 ) and f t (0) /(n 0 I 2 /rT 2 b 3 ) vs. 5 for 
0=5. In this graph, the forces are normalized to the magnitude of the total 
force exerted by the isolated wire on a test current at the corresponding 
point. Hie relative smallness of the tangential force is obvious. 

In the small 0 (0 « 1) regime, we find from Eqs. (2.3.27-31t 


H_I 

f t (0) " *b 


“V P 2 r i / 1 2 y-i'\ n 2P l , , V?' 

f n (0) = ~T~ T l. TTR 1 ) + v T- 3 -J e - yH ln T' 

TTb 

A 


<L__ 8g t 

16 V 3 tv ufr • 


(2.3.34) 

(2.3.35) 


We notice the following features of these expressions: 

1. The forces are in the same directions as they were for large 0. 

2. For given fixed values of b and F, the forces are much smaller for a 
given small value of 0 (say 0 - 1 /5) than they were for a corresponding large 
value of 0 (say R = 5) . This is simply because in the small-0 case the force 
distributions extend into the melt. 

3. The logarithmic term in f indicates that the expression (2,3.34) may 
not be valid for large § . 

4. For moderately large 5’s, the tangential force may be relatively large, 
even larger than the normal force. 

5. For fixed 6 and very small b, the shaping force f only increases 

2 n max 

as 1/b; while it increased as 1/b in the region where b was larger(Eq. (2.3.32)). 
Thus it seems probable that there is a point such that moving the wires in 
further to get a larger force will not be worthwhile, since the problems 
arising from any mechanical instabilities will be great, and the gain in 
shaping forces not so large. These remarks sufficiently describe the small-0 
situation so graphs of the force distributions will not be necessary, particu- 
larly since the regime where Eqs. (2,3.34-35) hold is not of great practical 
interest. Values of fields and forces for intermediate values of 0 could 
readily be obtained numerical lv if needed. 





FIG. 2.3.6. NORMAL AND TANGENTIAL 
MAGNITUDE OF FORCE ON 
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:e densities, normalized to 

r CURRENT AT ' - 0, FOR B - 5. 
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2.4 Determination of Degree of Shaping 


2,4.1 Introduction 

In this section, which for greater clarity will be divided into 
several subsections, we will develop expressions for the electromagnetic 
work done in changing the ribbon thickness, for the surface energy change 
with thickness, and for the shaper current required to shape the ribbon to 
a selected thickness. Power required for shaping is discussed in the fol- 
lowing section. We shall consider first an isolated ribbon — a strip of 

liquid metal of fixed cross-sectional area extending a long way in both 

( 2 ) 

directions from the shaper — and then the Me Donne 11-Doug las RMGS setup , 

where the ribbon is pulled from a sphere of molten Si near the shaper. 

Aside from its analytical simplicity, the isolated ribbon case is important 
for several reasons: 

1. The critical thickness- limiting factor for the isolated ribbon is the 
sharpening of the melt curvature at the ribbon edges. It has been necessary, 
owing to the complexity of the problem, to ignore this in the RMGS calcu- 
lation^ ^ ; so the isolated ribbon provides an oppositi- limit approximation 
to the true intermediate situation. 

2. In terrestrial experiments^, gravity will tend to inhibit shape changes 
along the meniscus between the ribbon and the molten silicon supply. Thus, 
the cross-sectional area at the shaping level will be difficult to change 
and the shaping should be rather similar to that in the isolated ribbon. 

3. Suggestions have been made for shaping thin ribbon by passing it through 
a series of progressively narrower coils, or by preliminarily shaping it 

by some other method and then finally thinning it electromagnetic ally. The 
isolated ribbon is clearly a good approximation to these situations. 

The ultimate thickness of the ribbon is conveniently calculated by 
the energy balance method. As the ribbon becomes thinner, the change in 
surface free energy to shape it further increases — that is, it becomes more 
and more difficult to shape the thinner it gets. Meanwhile, the normal 
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electromagnetic force on the ribbon Is less, the thinner It Is, and there 
will be less and less electromagnetic energy available to do the work of 
further shaping. The thickness at which the free energy required for an 
infinitesimal further thinning is just equal to the electromagnetic work 
required to do this thinning will be the actual thickness of the ribbon In 
a given field. We first determine this ultimate thickness for a flat ribbon 
of finite width with circular-cylindrical terminations, such as discussed 
In previous work^*^. 

2.4.2 Surface Free Energy for Isolated Ribbon 

We consider a ribbon of width W and thickness t, with curved 
euus of semicircular cross-section of radius t/2, as shown in Fig. 2.4.1. 

Most of the increase In surface energy on shaping comes from the end regions, 
which is why they must be included. If the central portion of the ribbon is 

Indeed flat, this configuration with the semicircular ends will be that of 

lowest free energy. We will assume that the shaping takes place uniformly 
over a region of height H In the pulling direction; the shape of the normal 
force vs. distance curve along this direction (Fig. 2.3.4) ensures that this 
approximation (made Implicitly by Gaule / and Pastore^^) will be reasonably 
good. Then If we consider the ribbon being shaped from a width W^ and thick- 
ness t ^ to a width W^ and thickness the change In surface free energy is 


A E 3 - v H[(2W 2 + nt 2 ) - (2W X + rrt^ ] 
with the constraint 

2 


“l*! + 4 C 1 


W 2 t 2 + 4 C 2 


(2.4.1) 


w'u^re v is the surface tension. If we let t 2 ~ t^ - ^t ( 80 At is a positive 

quantity), and take the limit of AE /At as At goes to zero, we find 

s 

2vHW 
t 


dE 
s 

dt 


(2.4.2) 
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FIG. 2.4.1. ASSUMED CONFIGURATION AND NOMENCLATURE FOR DETERMINATION 
OF DEGREE OF SHAPING OF MOLTEN RIBBON BY HAIRPIN COIL 
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2,4,3 Electromagnetic Pressure 

Although the eddy current force 1« a volume force, other dts~ 

(1 2 11 14) 

cussions * 9 ’ have assumed it could be considered as a pressure at the 

surface of the melt. For purposes of comparison with the volume force 
equations which follow, we will discuss this approximation briefly. If p(t) 
is the electromagnetic pressure exerted when the ribbon thickness Is t, the 
electromagnetic force on the shaping region Is p(t)WH, neglecting the small 
total force on the semicircular regions. Then the work done on one half of 
the melt In compressing It from thickness t/2 to (t - At)/2 Is |5WHAt/2, where 
p is the average pressure exerted over the range At/2. (Also W Is the 
average width, but we may ignore the weak dependence of W on t) . If At is 
small, p is substantially constant between t and t - At, and the total electro- 
magnetic work done on both faces is 


AWk * pWHAt 
or in the limit 


dWk 

dt 


pWH 


If ve imagine the ribbon being squeezed thinner and thinner by 
the field, the process will continue until the work required for an 
infinitesimal further reduction in thickness is just equal to that available, 
that is when 


dE s dWk 
dt + dt 


If this occurs at a thickness of, say, t f , then shaping stops when 

P(t f ) = t L 


.( 1 ) 


(2.4.3) 


(2.4.4) 
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2,4,4 Electromagnetic Work -Volume Forces 

If we describe the situation in terms of f^, the normal component 
of the volume force on the melt, then the pressure in the above equations 
is replaced by the integral of — in particular, 

■ WH f f n (y)dy (both Laces) . (2,4,5) 

'—CD 

This formula was derived by Gaule' and Pastore^ 1 ^ ; we obtained it independently, 
but omit the proof. The lower limit on the integral should more properly be 
-t/2, but because of the rapid falloff of the force inside the melt, taking 
the limit to ~<x> is an excellent approximation as long as the ribbon is at 
least 3 or 4 skin-depths thick. Ag-<n, the small contribution from the semi- 
circular regions has been neglected. 


2.4.5 Isolated Ribbon — Large B Approximation 


To proceed, the integral in Eq. (2.4,5) must be evaluated. This 
can be readily done in the regime where the parameter B> the ratio of the 
distance of the wire from the melt to the skin depth, is large. The type of 
approximations required have been already discussed (Section 2.3); we only 
have to generalize slightly to obtain approximations valid for y not equal 
to zero. On the assumption that B is large enough that only the first terms 
of the asymptotic expansions contribute, we get from Eqs. (2,3.11) and (2.3.12), 
after expanding the square roots, 


aud 


E z (x,y) 

H x (x,y) 


-I i+i 0TKl*i) 

TTFT- e 


-I _1 8T]( 1 1 i) 

~ nb 1 +p 2 


(2.4.6) 


(2.4.7) 


k 


* 
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where we recall 5 » x/b and 7] ■ y/b. Higher-order terms can be obtained If 
desired in terms of error functions of complex argument, but a more .^.factory 
procedure for obtaining the deyee of shaping for moderate 0's is dascribed 
in the next subsection. Then since 




we find 


f n * T Re (E z H x > 


f n ( 5"°> y) “ Wfi em 


or 


f n (?-0,y) - f nQ e" 2( l y/4 l ) 


(2.4.8a) 


(2.4.8b) 


in an obvious notation. In agreement with our previous approximation we only 
need f n for % « 0, so for the integral of gq. (2.4.5) we obtain 
,0 

f n d y - *n0 6/2 

_00 


Thus shaping may be considered to proceed until 


2vHW _ WHf n0 6 



(2.4.9) 


(2.4.10) 


with the constraint t f + 2b f « 2B, where 2B is the distance across the shap- 
ing coll (Fig. 2 .4.1). Eq. (2.4.10) may be rewritten 


! K b . 


K(B - 


: -± -) 2 


(2.4.11a) 


where the parameter 



(2.4.11b) 
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has the dimensions of an inverse length. Inserting the value v ' appropriate 
to liquid silicon, y * 0.72 N/m, we find 


.. 4.52 x 10 7 -1 

K » m 

1 

where I is the peak coil current in amperes. Sinca t f will be a fraction of 
a millimeter and since the rms current will commonly be given, we also write 


2,26 x 10 H 

(I V 

rms 


-1 

mm z. ^ 


K - — -y? 

where K, ■ 150 A/ (ram) ^2, 

1 

Eq. ( 2 . 4.1 la) is readily solved for 


— , say, 

v rms 


t 


f 


- 2B 



BK - A [ 2BK 'i 
BK J 


(2.4.12) 


Alternatively, we can consider we know the ribbon Lhickness we want to obtain 
and solve for the current lequiied to achieve it, given the size of the 


shaping coil: 

I 

rms 


B- t f /2 L 

K l"~tjr72 " K i T^2 


(2.4.13) 


where the lengths B, b, and are all to be expressed in millimeters . We 
repeat the conditions under which these equations can be used to predict the 
degree of shaping or current required: 

1. g * b/6 should be large (say 5 or more). 

2, 6 must be considerably less than t (say no mure than t/3) . 

These requirements may be met by operating at a high enough frequency; ve shall 
see in the following section that total power required increases gradually 
with frequency, but even at the lowest frequencies liable to be used, the 
approximat ions are likely to be satisfactory in any practical situation. It 
should be noticed that in this regime the current (though not the power) required 
for shaping to a given thickness is independent of the skin depth, and thus of 
the operating frequency. 
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To provide some numerical examples, if one wants a ribbon 0.4 mm 
thick, then If b is about 2.8 mm, as in the smallest coll discussed in Reference 
1, the rms current required is around 660 A, while if the wires can be brought 
in so that b is only .25 mm, the current required is 59 A. Clearly the trade- 
off of greatest importance is that between getting the wire close enough to 
get the current down and having it large enough to carry the current without 
melting. (Note that b is the distance from the melt to the center of the wire). 
Of course if the wire comes extremely close to the melt, the presence of the 
eddy currents will affect the current distribution in the wire; we will not 
pursue this extremely complicated problem, but just point out that Lens's Law 
shows that this effect will reduce the force on the melt, and that a somewhat 
higher current will be required for a given ribbon thickness than would be 

expected on the basis of the theory presented here. A sharp-pointed con- 
( 14 > 

centrator may be a good vay to get tne current close to the mel~; it will 
sti 11 clearly have to reach at least as close as the distance b, though. 

The work of Gaule' and Pastore^^ has been used^to estimate 
power and current requirements for electromagnetic shaping. Their approach 
and assumed physical situation are rather different from ours, but a comparison 
Insofar as is possible shows the following: 

1. They do not discuss the distance of the coll from the melt, and because 

of the nature of their approximations apparently did not realize the importance 
of this parameter. 

21 

2. They assume H (0,0) = — , and 6 * t/4. From Eq. (2.4.7), we see 

x c 

H (0,0) = — . By comparing these equations, we see that their assumptions 
x nb 

amount to taking 0 * 2/n. But otherwise their equations are similar to those 
based on the asymptotic expansions discussed above. Thus it appears that they 
are using asymptotic formulas in a region (0 not »1) where they do not apply. 

3. The fact that they considered Ge, with a surface tension about 0.6 N/m, 
rather than Si; a factor of 2 error in their equation for the pressure; and 
their other approximations enabled Gaule' and Pas tore to find the rather small 
value of 12 A required for shaping to .3 mn thickness. 
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2.4,6 Isolated Ribbon — Arbitrary B 


By evaluating an integral numerically, we may extent the results 
of the previous section to smaller p's and determine the range of validity 
of the asymptotic approximation used there. 

Eq. (2.4.5) may be rewritten 

dWk WH n n a r 0 

— = — — Re j E z (0,y)H x *(0,y)dy 

_Q0 


(2.4.5') 


Introducing the integral expressions (eqs. (2.3.11) and (2.3.12)) for 

E and H , we may write this as 
z x 


where 


dWk WHoo) I 


2 0 


, V > Im dy I I 
dt 2 V TT J 


It * 1 dk e 
1 *0 


-kb *r /Y.2 +■ at 2 y 


k * r A 2 


(2.4.14a) 


(2.4.14b) 


and 


I, . ■' Jo «- qb - ^ 


(2.4.14c) 


q i /q 2 - 0f2 

Changing variables of this triple integral to H * y/b, u * 6k and v = ftq, and 
doing the T| integral, which is a simple exponential, we find 


dWk __ -HWu 0 I F(g) 
dt 4 17^6 2 


(2.4.15a) 


where 


CD CO 

r 


F(0) du dv e " e(u+v ' Im f(u,v) 

0 1 0 


(2.4.15b) 


and 


f(u.v) = 


-4/yZ -21 


(u + /u2 + 21) (v + Jv l - 21) (/u* +• 21 + /v2 - 21) 


(2.4.15c) 
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The quantity f(u,v) can, if desired, be separated algebraically into its real 
and imaginary parts, but it is easier and faster computationally to let the 
computer take care of this. The function Im f(u,v) varies smoothly and is 
generally well-behaved in the first quadrant of the u,v-plane. It has a 
value of unity at the origin and falls smoothly with finite slope in all 
directions from there. The function becomes slightly negative near the u-axie 
for u greater than about 0.75, but this introduces no difficulty in evaluation 
of the integral. It is simpler for comparison with previous calculations to 
change variables again to s ■ 0u and t * 0v; then 


“ . -“V 2 0(B) 

dt 4 TT^b* 


(2.4.16a) 


G(0) - ! f da dt e~ (s + ^Im f(a/B,t/e) , (2.4.16b) 

v 0 “0 


with f as given above. Comparison with Eqs. (2.4.5) and (2.4.9) shows our 
previous asymptotic approximation amounts to taking G(0) * 1, independent of 
&. 

Since the integrand of (2.4.16b) is in the form of a double Laplace 

transform, it seems natural to try to evaluate by means of a product Gauss- 

( 21 ) 

Laguerre numerical integration routine . Such a routine was written and 
used to evaluate G(8) in the range 0.5 £ 0 £ 44. Fig. 2.4.2 shows results 
for 0 between 0.5 and 6. At 0 = 44, G(0) =* 0.9770. The evaluation used the 
first 16 x 16 points of a 28 x 28 Gauss-Laguerre scheme. (Because of the 
generally decreasing nature of the integrand, points further from the origin 
do not contribute appreciably to the integral, and it is a waste of time to 
evaluate them) . The accuracy of the integration routine was tested by trying 
some integrands which could be integrated numerically, and also by redoing 
the integral of interest here for several 0’s by 16 x 16 point Gauss-Legendre 
quadrature in 36 or 81 subintervals. In the latter procedure, appropriate 
upper-limit cutoffs to tue integral must be Introduced; there is no problem 
In doing this so that the remainder is negligible. It was found that the 
simple Gauss-Laguerre procedure could be relied on to 5 decimal places for 
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0 Ss 2.5, and gave an error in G(6) of only 0.77 percent for 6 - 0.5. The 
error as a function of A is graphed in Fig. 2.4.3. The value returned by the 
Gauss-Laguerre quadrature program is invariably slightly higher than the true 
value. The quadrature scheme will be rather Inaccurate for very small B's, 
but the small 8 region is not of much physical Interest. 

It is apparent that the asymptotic approximation is not very accurate 
for 8*8 in the range discussed earlier. For the problem of shaping at con- 
stant volume discussed there, all that needs to be done to make the theorv 

1/2 

correct for any 8 is to replace b f in Eq. (2.4.10) by b f /(G(g)) . Thus, the 

r * - 1/2 
current required to achieve a given degree of shaping will be a factor (G(B)) 

larger than would be calculated from the asymptotic expansion. For 8 = 2/tt, 

the current required for shaping to a given thickness is 2.1 times larger than 

the asymptotic value, 

2.4.7 Degree of Shaping in RMGS Configuration 

( 2 ) 

In a recent study at McDonnell-Douglas of potentialities for 
space manufacture of Si, a technique termed Ribbon from Melt Growth in Space, 
or RMGS, has been described. It is discussed briefly elsewhere in the report; 
the salient idea as far as the present discussion is concerned is that the 
ribbon is drawn through the shaper from a freely floating sphere of molten 
Si, The constraint of constant cross-section then no longer operates; one 
may look at the process as one in which the excess Si is squeezed back into 
the sphere when the coil current is turned up to thin the ribbon. Since the 
sphere radius is large, one might anticipate that the free energy change on 
shaping would be less than that calculated previously. We wish to assess the 
effect of this change in configuration on the degree of shaping for a given 
coil current. 

An analysis has recently been made by R, Rochat of McDonnell- 
Douglas of the expected shape of the meniscus from the ribbon to the sphere 
and of shaping pressure, assumed applied uniformly over a height H, required 
to shape the ribbon to a given final thickness. Numerical methods were 
required for this analysis and the results have been presented graphically. 
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FIG. 2.4.3. rRROR IN CALCULATION OF G(P) BY 16 x 16 ORDER 
PRODUCT GAUSS- LA GUERRE QUADRATURE 
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We have not had an opportunity to check this work independently, but find no 
reason to doubt Its correctness. The effects of the curvature at the edge 
of the ribbon have not been included, but this may not be a severe approximation 
in these circumstances. One parameter entering the calculation is the ratio 
of the height of the region being shaped to the total vertical height of the 
meniscus region. The present discussion will be most appropriate when this 
ratio is not too large — say 1/2 or less; however, the results of the surface 
stress analysis do not depend strongly on this parameter. 

The principal result of the analysis is that the shaping pressure 
is given, instead of by Eq. (2.4.4), by 

p = yf/H , (2.4.17) 

where f is a dimensionless function* of the ribbon thickness t, the sphere 
radius R g , the ratio N discussed above, and H. The principal properties 
of the function f for our purposes are 

1. f is always greater than one. 

2. for any reasonable set of parameters, the value of /7~, which we 

will see to be the parameter of greatest Interest, will lie between 1 and 2. 
Generally it will be rather close to 1. 

From Eqs. (2.4.3), 2.4.16) and (2.4.17), we find the peak current required 
for shaping to thickness t under these circumstances is given by 

I RMGS = V p, 0 HG(3)y ' (2.4.18) 

It is interesting to compare this to the current required for shaping of an 
isolated ribbon, I R , which can be found from Eqs. (2.4.2), (2.4.3) and (2.4.16) . 

We obtain 


I c 1/2 

RMGS _ _X_ 

T “ V 2H G ^ 

R 


(2. -.19) 


Elsewhere in the report, f denotes frequency. This dimensionless 
function f will only be used in this subsection. 
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f/G will be greater than unity. t/2H might be small, but if H is large so 
are b and, consequently, I D . To pursue this further, we need a specific 
relationship between H, the height effectively shaped, and b, the distance 
of the wire from the melt. A reasonable approximation may be obtained by 
averaging the normal force, given by Eq. (2.3.32) for large 0's, over §, 
the normalized distance along the ribbon. One obtains 


where 


H = - 2 b f (8) 


r(8) 



(2.4.20a) 

(2.4.20b) 


Now let us specialize to the case we have frequently discussed, 
where 6 = t/4. (One can readily generalize to other ratios; tiis is 
discussed qualitatively below). Then we can express the current ratio 
almost completely (except for *Jf ) in terms of the dimensionless parameter 

9: 


I 


RMGS 


$(0) /f 


(2.4.21a) 


where 


•(B) =(4/n8rG) 1/2 


(2.4.21b) 


$(0) is plotted in Fig. 2.4.4, We see that for small enough 0--0 less than 

2.5 for as l--the ratio L u . c /I_ will be greater than unity; so for small 

KMGo K 

0 it will be more difficult to shape the ribbon when it is attached to the 

spherical melt. For larger 0‘s, though, there will be a diminution of current 

required for shaping to a given thickness from a given distance when the 

presence of the spherical melt is considered. 

1 will also depend on 0 of course; so more current will be 
R 

required when b is larger. Again assuming t/6 * 4, we can write 
1 R = tt0 (26’v/uq) 1/2 ; 


so 


RMGS 


= TT 


(26y/li 0 ) 


1/2 


0 §(B) 


(2.4.22) 
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0 $ (0) is also plotted in Fig. 2.4.4 and is seen to increase slowly with 0 . 

Thus the advantage to be gained by bringing the wire in close to the melt 
is not extremely large; on the other hand, the penalty for moving the wire 
further away is not so great. To give a specific numerical example, suppose 
6 = 100 M-m. Then the rms current required for shaping a 400 ribbon 
is given by 

X RMGS “ 23,8 9 t(0) /f A • (2.4.23) 

Taking N /f - 1.03 as a not untypical value, we find the current required 
for various distances from the shaper to the melt as in the following 
table : 


TABLE II. rms CURRENT AT 21 MHz REQUIRED TO 
SHAPE Si RIBBON TO 400 u m THICKNESS 


Sbaper-to melt 
distance, b(^m) 

200 

500 

1000 

1200 


rms current 
(A) 

59 

73 

95 

102 


If we had assumed a ratio t /6 - t, say, instead of 4, we would obtain a 
1/2 

factor of (4 /t) on the right side of (2.4.23). Thus some saving in cur- 
rent required can be obtained by reducing 6 . But the only way to do this 
is to increase the operating frequency. Since 6 ^ 0 ) I 01 uu so the 

frequency must be increased by a factor of 16, to over 330 MHz in the 
situation described above, to cut the current required by half* The ac 
resistance of both the shaper coil and the melt will also increase with uu; 
the combined effects of all these changes are discussed in Section 2.5. 


* or less, if the frequency dependence of 6 $ (0) is considered. 
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2,4.8 Degree of Flatness 

We have assumed to this point that the ribbon was essentially flat, 
its cross section at the plane of the shaper being a rectangle terminated by 
two semicircles as in Fig, 2.4.1. While the physical situation surely 
indicates that the ribbon will be fairly flat, it would be useful to have 
some measure of what the departure from flatness might be. In this section, 
we estimate this quantity for the case of large 3 where the departure from 
flatness is most likely to be significant. To make such an estimate, we have 
to make some assumption as to the shape of the ribbon cross section, since 
we do not know what the shape of lowest total energy will be. We will take 
the cross section to be an ellipse of the same area as the flat ribbon. Then 
by trying ellipses of different eccentricities, we can determine the elliptical 
shape of least total energy, which should be good approximation to the actual 
shape taken up by the ribbon at the plane of the shaper. 

The geometry of the situation is illustrated in Fig. 2.4.5. The 
equation of the melt surface is 



where the area of the cross section is 

2 

A * tt pq ■ Wt + TTt /4 

Since A is fixed, if we select a trial value for the semi-minor axis p 
(presumably a little more than t/2> , q is determined; so there is just one 
free parameter to vary. 

Again assuming the shaping takes place uniformly over a height H, 
the change in surface energy on going from the flat to the elliptical cross 
section is 


AE = 


"y H 1 4q E(/ q 2 - p 2/q) - (2W • nt)1 


(2.4.24) 
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FIG. 2.4.5* DEFORMATION OF FLAT RIBBON OF WIDTH W, THICKNESS 
t, INTO RIBBON OF ELLIPTICAL CROSS SECTION WITH 
SEMIAXES p and q 


i 



47 


i 


i 


( 22 ) 

where E denotes a complete elliptic integral . The corresponding change 
in electromagnetic energy may be approximated by assuming that each small 
element A* of the melt is part of a very wide flat plate of thickness t(r) 
and distance b(t) from the shaper. Since the shape of the melt only ctu.^es 
gradually, this should he a good approximation. Then the change in electro- 
magnetic energy on going from the flat to the elliptical configuration is 
obtained by calculating the energy change for each region of width Az 
using Eqs. (2.4.5) and (2.4.8), and Bumning over all the strips from z » -W/2 
to z * +W/2. For consistency with the previous approximations we again 
neglect, if possible, the electromagnetic forces on the very edges of the 
ribbon. Thus we have 

AWk - 2H [w 0 I 2/ An 2 ]/o iz ^2 dt< fb(t ' )l2 • (2 - U 25 > 


where 


u = min [W/2, q] 

The current I is fixed at the value that vas calculated necessary to shape 
a flat ribbon of thickness t, so it is given by Eq. (2.4.10): 

PLqI 2 /4tt 2 « 2yb^ 2 /t , 

where b^ denotes the value of b for a flat ribbon. What we want to do is to 
calculate the total energy change A E total * A^ 8 *■ AWk as a function of the 
amount of "bulge 11 in the ribbon (l in Fig. 2.4.5), and the other parameters, 
and to see whether there is some value of t for which the energy change is a 
minimum and negative — if there is, a configuration similar to that elliptical 
one will probably be taken up by the ribbon The integrals in Eq. (2.4.25) 
can be done analytically, but since the results are rather complicated, we 
shall Just present some results graphically. It will be convenienc to divide 
both AE g and AWk by the factor (2 y WH/tt) , and to compare the dimensionless 
quantities resulting. If we also take as a dimensionless independent variable 
the fractional bulge at the ribbon center* or ■ 2l/t, we find there are only 


* This a should not be confused with the at that vas used in section 2.3 
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two dimensionless parameters In the problem: 

1. s * rrt/4W, a quantitv depending on the aspect of the flat ribbon Both 
AE 8 and AWk depend only very weakly on this parameter, but we will assume a 
large but finite aspect ratio, W/t * 150, for purposes of calculation. 

2. r * 2b^/t, a measure of the relative size of the shaping coil. AE S of 
course is independent of this parametar. For our usual assumed conditions, 
where t « 46, r is related to our usual parameter 8 by r « 8/2. 

Some calculated results are shown in Fig. 2.4.6. It appears that 
although there is no small change from flat ribbon for which the departure 
results in a decreased total energy, this is primarily because a large bulge 
is necessary for the surface energy to decrease, and it does appear that for 
some increase in thickness in excess of 30 percent, an elllotical cross 
section will be more stable. We have not made calculations to determine the 
minimum energy configuration for this region though, since with further 
increases in thickness, we are getting rapdily into the region where the 
asymptotic approximation breaks down. Even *f more elaborate calculations 
(which could be carried out numerically) showed that the total energy change 
was never quite negative, but only reached a minimum, the presence of such 
a low-energy configuration could be important in determining the fitaoility 
of the system against vibrations or other mechanical perturbations. Finally, 
we wish to re-emphasize that there are liable to be changes in ribbon shape 
between the shaper plane and the level of the freeze, which we have not con- 
sidered in this section. 


2.5 Power Requirements 

In this section we estimate the cotal power needed for shaping to 
a specified thickness, with a coil >f given dimensions, in both the isolated 
ribbon and RMGS situations. 
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a= 2 L/ t 


FU,. 2.4.6. CHANGE IN ELECTROMAGNETIC, SURFACE, AND TOTAL 

ENERGY WITH RELATIVE ELLIPTICAL DISTORTION OF MELT 
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The total steady power needed at the input to the shaper coil 
will be the sum of the power dissipated in the melt, P m , and that dissipated 
in the coil, P c . The power dissipated in the melt is given by the real part 
of the integral of the Poynting vector over the ribbon surface, Eq. (2.2.13) 
In the asymptotic (large 8) regime, we find, using just the first terms of 
(2.3.20) and (2.3.21), 


P m = 


2 -1 
» I fWB 
^o rms 


(2.5.T 


where we recall W is the ribbon width and f is the operating frequency. Use 

of the asymptotic values for the fields is not as severe an approximation as 

it was for determining degree of shaping, and (2.5.1) will be correct to 

within 10 percent at 8 = 5 and more accurate for larger B's. 

If we take the coil to have a circular cross section of radius a, 

(23) 

then the power dissipated in the coil is given by 


I R 
rms ac 


where 

R ac = (a/26 c )R dc 


(2.5.2a) 


(2.5.2b) 


*dc - 


L/rra^) 


(2.5 2c) 


where p c is the resistivity of the coil material and L is the total length 
of the coil. 6 C is the skin depth in the coil, determined by and the 
operating frequency, and we have made the generally very good assumption that 


a >2 


c * 


The power dissipated in the coil will frequently be comparable to 


that dissipated in Che melt so ic must be included in any realistic estimate. 
Adding (2.5.1) and (2.5.2), we obtain 


^O^rms 


f 6 


V 4 L_ *£- 
b 2a 6 _ 


(2.5.3) 


Inserting the value of the current obtained from Eq. (2,4.10), we find for 
the isolated ribbon 


IR 


An^vffib^t 


W 

b 


L 6c“ 

2a 6 


(2.5.4) 
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1/2 - 1/2 

In this form it is clear that * f , or « 6 . Thus for a specified 

ribbon thickness and coil dimension, the power required may be reduced 
somewhat by running at a larger skin depth. Of course with too large a 
skin depth the nornu.1 shaping force on one face of the ribbon will be re- 
duced by the force coming through from the other side, and the shaping 
effect will be reduced. Presumably the optimum 6 for given conditions could 
be established by a numerical procedure such as that described in Section 2.6. 
Furthermore, we can see from the discussion of Section 2.4.6 that the current 
required will increase somewhat with decrease in frequency. Selection of 
suitable operating voltage and frequency cannot therefore be made quite 
independently. 

The necessity of avoiding contact between the ribbon and the 

shaper will doubtless lead to some minimum clearance that can be allowed 

between the shaper and the melt. While this will naturally depend on the 

other parameters of the problem, let us assume that it can be taken as an 

independent design variable and call it c. Thus c = b-a. In the isolated 

ribbon situation, it is easy to show that there is an optimum value of the 

ratio c/b for which the power required is minimum. If we assume, as in 

some of the McDonnell- Douglas calculations, that L * 3W, and further assume 

-8 

that the coil is copper, of resistivity 1.7 x 10 ohm m, then the optimum 
value of c/b is 0.704, and roughly 407« of the power is dissipated in the 
coil and 60% in the melt. If we assume our usual operating conditions 
f * 21.1 MHz, t * 400 p.m, W = 6 cm, then the total power required is simply 
given by 


P t IR (kW) = 22.06 c (2 5.5) 

where the clearance c is in millimeters. For any reasonable value of c, 
this is quite a substantial power requirement, enough to make electro- 
magnetic shaping of an isolated ribbon relatively unattractive. 

In the RMGS situation, if we take as unity the function f dis- 
r ised in Section 2.4.7, we have 

* P c IR (i(P)] 2 , (2.5.6) 


where $(0) was given by (2.4.21b) and we must remember that we have already 
specialized to the case 6 = t/4. In this case we find that rather than 
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there being an optimum cA>, for fixed c the power required approaches 
a minimum as b-* 00 . Explicitly, 

P^^Onin) = 4nYfWJ(l+(L/2W)(6 c /6)) . (2.5.7) 

For the same conditions as used above for the isolated ribbon, we find 

RMGS 

that this minimum power is only (min) » 1.4 kW. Of course the 

actual power required will be more than this, since b cannot be too large 
for the current to be reasonable-- a few values are given in Table III--; 
but in any event the total power requirement in these circumstances seems 
quite moderate, and much less than for the isolated ribbon in similar 
circumstances. Thus although the current required in the RMGS configura- 
tion is higher than that which might be expected from naive considerations, 
the power requirements are such that the process appears from this stand- 
point rather attractive. Of course the interaction with the meniscus 
region, which would certainly increase the power required, ha,? not been 
included. A realistic calculation of this effect appears extraordinarily 
difficult , though. 


TABLE III. POWER (IN kW) REQUIRED FOR SHAPING RMGS 

CONFIGURATION AS FUNCTION OF DISTANCE b FROM 
MELT SURFACE TO CENTER OF SHAPER AND OF CLEAR- 
ANCE c BETWEEN SHAPER AND MELT. 



0.25 

0.5 

1.0 

1.5 

2.0 

.8 

1.84 

2.20 




1.2 

1.66 

1.79 

2.99 



1.6 

1.59 

1.66 

1.99 

5.61 


2.0 

1.54 

1.59 

1.77 

2.30 


2.4 

1.51 

1.55 

1.67 

1.92 

2.80 

2.8 

1.50 

1.53 

1.62 

1.77 

2.12 
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2.6 Effects of Freezing Interface 
and of Finite Ribbon Thickness 


It will be shown in Section 2.7 that the liquid-solid interface 
will typically be veil away from the shaping region; thus treating this region 
as all liquid silicon will generally be a very good approximation. It is 
possible, though, that the Interface might be brought closer to the shaping 
level by some auxiliary cooling method, and it would be of Interest to know 
quantitatively what changes the presence of the solid region would make to 
the fields, forces, degree of shaping, etc. In this section, we outline a 
numerical procedure which would enable one to determine the electric field 
distribution in the situation described. From this, the magnetic field, 
forces, etc., can be estimated by standard numerical techniques. I there 
is negligible time delay in the shaping coil, a strict boundary condition at 
the center plane of the ribbon (solid or liquid) Is that the electric field 
be zero. By adopting this boundary condition, we are also enabled to investigate 
effects resulting from the ribbon being thinner than a few skin depths (though 
at least one skin depth thick) by such a numerical technique. While computer 
subroutines have been written to handle most aspects of the numerical procedures 
involved, lack of time at the conclusion of the program has prevented us from 
running these programs to any great extent. 

(24) 

The method adopted is a modification of one due to Stoll. This 

method is essentially one of finite-difference approximation to the scalar 
Helmholtz equations for the complex electric field. Stoll has ingenuously 
modified this equation, however, to take the magnetic field boundary conditions 
into account, by assuming a smoothed spatial variation of the magnetic 
permeability. The method is applicable to fields which depend on two dimensions 
only and to regions with rectangular boundaries. Extension to more general 
regions might be possible, but would be quite difficult. On the other hand, 
the method is readily applied to any number of subregions of different materials 
and to any number of current sources, either distributed or infinitesimal. 

To solve the Helmholtz equations by a finite-difference method, it 

is necessary to fix the field on some suitable boundary of the region within which 

(24) 

a mesh of sample points is to be defined. In the problems discussed by Stoll, 
it is possible to choose a boundary on which the electric field is negligibly 
small, and may be taken to be zero. Even in these cases, many points outside 
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the eddy current region must be included in the mesh, though. In ou situ- 
ation, where there is a current source in the vacuum region outside the melt, 
including enough points in the vacuum region that the field be negligible 
is out of the question. We can, however, find a boundary in the vacuum 
region where the real part of the field is negligible and the imaginary part 
is close to that of the problem previously solved with the ribbon assumed 
thick and entirely liquid. Thus by having solved the previous problem, we 
can obtain a suitable starting point from which we can bootstrap our way to 
a solution in the present more complex case. 

It is possible again, as in Sec. 2.3, to normalize all distances 
to b. the distance of the shaping wire fro© the melt. Thus the only param- 
eters entering the problem are the ratio of the shaper distance b to the 
skin depth in the liquid 6; - x Q /b, the relative distance to the solid- 

liquid interface from the shaping plane; P * P/P 8 > th® ratio of the resistivi- 
ties in the liquid and the solid, and the exciting current I. We will see 
that the numerical procedure will only be reliable, with a mesh of reasonable 
size, for rather large g*8 — say 6 ^ 10. The setup of the region over which 
the finite-difference mesh is to be described is shown in Fig. 2.6.1. We 

recall that £ » x/b and T) - v/b. By consideration of our previous results 

(24) 

and of numerical results on simpler problems, we judge that a mesh size 
of A? = *2 and ATj » ,2g - ^ in the liquid and solid, AT| * .2 in the vacuum 
out to T| a* 4 } and AT] * 1 for T] > 4 should be sufficient to determine the 
fields to within a few percent. A portion of this mesh is also shown in 
Fig. 2.6.1. 

From Eqs. (2.3.3c) and (2.3.10a), we can show, using the identity^ 


-v 


2 . 1/2 ... 

+ or ) /(k 


<k‘ 


a 2 ) 1/2 ) 


-1 


(k + (k* + 


«V /2 > 


( 2 . 6 . 1 ) 


that tu 2 unperturbed electric £ield in the vacuum region is given, for large 
8, by 


-U 0 *l 


2 2 

+ (Tj+1) 

1/2 In — s r. + x (u,v) ; 

? 2 4- (T|-l) 2 J 


2rr 


( 2 . 6 . 2 ) 
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FIG. 2.6.1. FINITE- DIFFERENCE GRID IN NORMALIZED COORDINATES 
FOR NUMERICAL SOLUTION OF A TWO-DIMENSIONAL EDDY' 
CURRENT PROBLEM 


where 


U = B(T|+1) , v = (? / (T|+l) ) 2 , and 


X (u,v) 



i 1-v 1+1 l-3v 

u (l+ v 2)2 " 2u ' TT+v2) 3 


By evaluating this expression for large 0's and by numerical calculation of 
an integral for x in a way similar to that of Sec. 2.4.6, we found, for 
instance, that for 0 £ 10: 1) the real part of x would be no more than 7% 

of the magnitude of the field at the origin at Tj = 14, 2) the imaginary 
part of x would be more than one order of magnitude smaller than the real 
part, and 3) Re x is no more than 5 % of the logarithmic term. The situation 
would of course be improved if we took the skin depth as that of the solid. 

Thus if we take T) = 14 as the right-hand boundary of the sample region, and 
fix the electric field there as that due to the logarithmic term (which will 
be the same regardless of P Q ) , we might anticipate reasonably accurate results 
(though possibly 5 to 10% off) in the region of interest near the ribbon surface. 
A similar analysis along the ribbon-vacuum interface (where the logarithmic 
term is zero) indicates that even for = 0, similar accuracy couM be obtained 
by going from ^ = 4 to ? = +6, again fixing E along these lines as the value 
of the logarithmic term. Interior to the ribbon, it will be satisfactory to 
take E » 0 on the horizontal boundaries as well as at the ribbon mid- plane. 

With the mesh described above, there are 1911 interior points at which the 
field must be determined if these limits on § and 7| are adopted. Of course 
it is prudent to make the limits variable in computer programing in order 
to allow for possible poor convergence or for future refinement, but the 
amount of computer time and memory required mount rapidly. 

By hypothesis the exciting current wire (which we shall take to be 
of infinitesimal size) lies on a mesh point. Since the mesh points are closely 
spaced, it will be satisfactory to assume that the interface plane P * 
passes through a row of mesh points. One parameter of the problem is eliminated 
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if we define all fields in terms of a reference field E D - pl/(2nfi ). Then 
the field at a given mesh point E Q is given in terms of fields at adjoining 
mesh points (Fig. 2.6.2), for each of the regions of Fig. 2.6.1, by 

the following equations: 
in region v ? , 

E 0 " (E 1 + E J + 25(E 2 +E 4 ))/52 ; (2.6.3) 


on the v 2 ~ v i boundary, 

E Q * (E 1 +5E 3 +15(E 2 +E 4 ))/36 ; 

in region v^ 

E 0 - (E 1 +E 2 +E 3 +E 4 -in6(?)6(Tl-l))/4 
in the liquid, 

B 2 (e 1 +e 3 )+e 2 +e 4 

£ - - - 5 

U 2 (l+r+ip z /25) 

in the solid, 

3 2 (e 1+ e 3 )^e 2 +e 4 

E 0 = 2 2 * 

2 (1+3 .g P/25) 

on the liquid-vacuum boundary, 

2&E 1 +2P 2 E 3 +(1+ME 2 +E 4 ) 

r _ _ - ; 

2 ((i+er+ir/25) 


(2.6.41 


(2.6.5) 


( 2 . 6 . 6 ) 


(2.6.7) 


( 2 . 6 . 8 ) 
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on the so lid -vacuum boundary, 

20E 1 +2e 2 E 3 +(l+B)(E 2 +E /t ) 

E 0 “ 2 2 

2 ((1+B) +KTP/25) 


(2.6.9) 


on Che solid-liquid boundary, 
e 2 (E +EJ+E +E, 

E 0 2 2 ( 2 . 6 . 10 ) 

2 (1+B )+l(l+P)B /25 


and at the point where solid, liquid and vacuum meet, 


2pEj+2S E 3 +(1+B)(E 2 +E 4 ) 
0 2 (l+8) 2 +i(l+P)P 2 /25 


( 2 . 6 . 11 ) 


together with the conditions at the points on the boundary: 


-i In 


S 2 ♦ (Vl ) 2 

5 2 +(T]-i) 2 


( 2 . 6 . 12 ) 


on the external boundary; and 


E = 0 


(2.6.13) 


on the boundary inside the riobon. In (2.6.5), the deltas are defined by 

6(x) = 1 if x = 0, 6(x) » 0 otherwise. 

There are a number of ways of looking at the problem of solving the 

very large set of coupled equations represented by (2.6.3) - (2.6.13). One 

might for instance, consider it simply as inversion of a very large sparse 

matrix, albeit a complex and not necessarily square one. A better approach 

seems to be to consider it a problem in adaptive control, or so-called "machine 

(24) 

learning". Such a process, described by Stoll, involves first of all the 

usual procedure of successive over- relaxation — that is, the mesh is repeatedly 
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scanned by columns and the value of E at each point modified according to the 
above equations, without regard to whether some neighboring values have been 
modified on the current scan or not. If this process converges, it can be carried 
on until some convergence criterion is satisfied; Stoll suggests the maximum 
difference in E's at any point between successive scans be less than some selected 
value . In the present problem, it might be better to have a more stringent 
condition applying to the region of greatest interest near the origin than 
elsewhere. The criterion might apply to the real and imaginary parts 
separately, or might Just be made to apply to the modulus. This simple 
process can be greatly speeded up by Increasing the field change after every 
scan by some (complex) convergence factor of modulus greater than unity. 

After every few s'-ins, this factor can be modified according to the results 
so far obtained. Implementing such an adaptive control technique for a com- 
plex field is rather tricky, but a surprisingly simple method which has worked 
successfully on a number of problems has been developed by Stoll. It is 
attractive to consider also the possibility of gradually varying some of the 
other parameters in the problem simultaneously. For Instance, one might 
gradually vary the resistivity ratio P from an initial value of one, for 
which the solutions have been determined in Section 2.3, down to its actual 
value of .05. We have prepared computer subroutines for evaluating the electric 
fields at the mesh points in the above equations, for cycling through the 
relaxation procedure with arbitrary convergence factor, for initialization, 
and for suitable input and output, but we have not programmed the adaptive 
control techniques just described. 

It might be possible to modify the techniques described in this 
section to estimate the fields in the RMGS configuration by replacing the 
meniscus by a right-angled bend. 



1 


r 




61 

2,7 Thermal Effects 

In this section we shall consider separately heating effects 
arising from the eddy currents in the silicon and those resulting from rhe 
current in the shaping coil. 

2.7.1 Eddy Current Heating in Si Ribbon 

Eddy currents strong enough to shape a thin ribbon will also be 
strong enough to heat it significantly. In this section we consider the 
effect of this heating on the position of the solid-liquid interlace and on 
the maximum steady-state temperature in the ribbon. 

The statement made in Section 1.3 that only steady-state temperature 
conditions need be considered requires some additional support in view or the 
presence of the liquid-solid interface, which cannot be moved vithou^ supplying 
the latent heat of 1.8 MJ/kg (Table I). From Table III, ve see that the 
power input will be around 2-3 kW. Anticipating some later results of this 
section, we find the amount of reradiated pover is at the most .5 kW, and 
might be somewhat smaller. The small dimensions of the ribbon, the high 
temperature of the melt, and the lower thermal conductivity of the solid will 
limit conduction losses (thic again may be justified a posteriori from 
results below); no most of 1.5 kW net iuput power should be available to 
move the interface. The initial velocity of the interface can be determin 
from 



? = L A u , (2.7.1) 

net V I 

-3 

where L. . , the latent heat per unit volume, is about 4.5 GJ m , and the 

^ -5 2 

ribbon cross-sectional area A is taken as 2.4 x 10 m . Thus u^ 1.5 

cm/s ec. Since this is much larger than anticipated pull rates, we 

anticipate that the interface, if initially at the shaper plane, would very 

rapidly move to its steady-state position. It is interesting to note that 

this po&ition does not depend on the latent heat, which only affects the 

dynamics of the interface. 
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Because the ribbon is so thin, the temperature variation across 
it can easily be neglected. Also since there will be little radiation from 
the ribbon edges, the temperature will vary but little across the breadth 
of the ribbon. Thus we have a problem of sceady-state heat transfer in one 
dimension-- the pulling direction, with radiation from the ribbon surface 
and heat produced within it. To avoid analytical complexity we shall 
ar* .-me that the heat is produced uniformly within a region of width H, as 
we assumed in discussing the degree of shaping. The bulk of the melt, \ here 
the temperature will ue fixed at T , the uelting temperature, will be taken 
as at a distance L from the plane of the shaper. The change in ribbon 
thickness near this point will be neglected; this can be shown to be insig' 
nif leant. At large distances along the solid ribbon, the temperature will 
approach the ambient, which since it will be relatively low we shall take 
as zero. The overall situation is thus as depicted in Fig. 2.7.1. 

The power radiated from a strip of ribbon of height Ax is gi<en 
4 

by P ra( j = 2W £e(T-T 0 ) , where £ is the Stefan- Boltzmann constant and T c 
the temperature of the surroundings. The heat transfer problem may be 
simply solved analytically ii we can make a linear approximation to the 
quadratic temperature dependence. When the temperature change from tUe 
ambient is small, the approximation 

(T-T 0 ) 4 « 4T 0 3 (T-r o ) 

is frequently used, but in the present situation, where the temperature 
change is large and T q ^ 0, this is clearly inadequate. The linear approxi- 
mation can still be used, though, ii we define a fictitious temperature T* 
such that 

4 *3 

T T , (2. 7 .2) 

★ 

where we have set T^ * 0, and where T is determined self- consistently in 

such a way that the radiation at temperatures ne&/ the actual maximum 

★ 

operating temperature is given accurately. If T is chosen such that 
radiation at the maximum temperature anywhere in the ribbon is given 
correctly then the overall radiation loss is slightly overestimated. The 
numerical prccedure for obtaining the maximum temperature will be seen to 
be quite straightforward. 


i 
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Thus in the four regions indictee! in Fig. 2.7.1, the heat 
conduction equations to be solved are simply 


d^T 

^ = c i T i * 0 * 

-L < x < -H/2 ; 


(2.7.3) 

= c i t ii + c 2 = 

0 , -H/2 < x < H/2 

> 

(2.7.4) 

d T m _ 0 

^r- = c i T m - 0 

, H/2 x < X Q ; 


(2 . " .5) 

b^ t 

IV 

C s IV 

» x 0 x v. « 


(2.7.6) 

= 2E e T* 3 /tK , 



(2.7.7) 

c s " 2E e s T * 3 / tK s * 



(2.7.8) 

c 2 = 0I 2 /(2nHb6tK) 

> 


(2.7.9a) 


where 


and 


or if Eq. (2.4.10) holds 

c 2 = 4TT 2 YfPH -1 K~ l («/t) 2 

For numerical work we will use (2.7.9b) with t * 46, f = 21.1 MHz, and 
H = nb/2. The implicit assumption of these equations that the interface 
position X 0 >H/2 will be seen to be well justified. 

Eqs. (2.7,3)- (2.7.6) together with the boundary conditions 


(2.7.9b) 


and 


T I = T m * 

X = -L ; 

(2.7.10) 

t = T ir a V* x “ 

dTn/dx, x = - H/2 ; 

(2.7.11) 

II " T IV* dT n /5x 

= dT m /dx, x = H/2 ; 

(2.7.12) 

III ■ T IV K ,T IlI /d ’ 1 ■ K . V*’’ * ' X o • 

(2.7.13) 

IV -°, 

1 

(2.7.14) 
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are readily solved. We shall not present the general solution, but Just 
give the temperature at the solid- liquid interface 

T m + 2R sinh sinh Z L 

T i “ T(x “ V * S sinh ( F + ZJ + cosh ( F + ^) 
and the maximum temperature in the ribbon 


T max * T ( x s = R (1 * cos ^ Z j^ + Tf (S sinh Z q + cosh Z Q ) . (2.7 

In these equations, 

R = TTYf/(4S eT* 3 ) , (2.7 

1/2 

S - (K e /Ke) = 0.63 , (2.7 

s s 


h ~ c i " L 


Z H SC 1 


1/2 

1/2 


H/2 


and 


1/2 

Z ■ c ' X 
c 1 o 


(2.7 

(2.7 

(2.7 


* 3 5 

If T =10 K, then with the other assumptions as before, R = 4.21 x 10 K 

1 **2 -1 

and c^ ' = 0.92 cm . Because of the asymmetry of the problem, T^ y 

might actually be slightly away from x = 0, but we can safely ignore this. 

To solve (2.7.15) and (2.7.16), we assume a value of T , calculate 

T = T /4 i/J , dete^aine R and 3, determine Z numerically from (2.7.15) 
max o 

assuming there is no supercool ing so T, « T , and calculate T from 

i m max 

(2.7.16). The true T will lie between that assumed and that calculated, 
v max 

and can be found with only a few guesses. We carried this procedure out for 
L = 4 cm and L = 0.4 cm (about the largest and the smallest values one might 
wish to consider ^^) and for b = 0.5 mm and 2 mm. Even for the smaller value 
of b, we found X « 0.3 cm for both values of L and T of 3700-3800 K, well 
above the silicon boiling point (at 1 atmosphere) of 3060 K. Thus it appears 
that the interface will generally be well out of the heating region, and 
also that some auxiliary cooling method will be necessary to keep the silicon 
at a reasonable temperature. Such cooling would bring the interface closer 
to, but probably not into, the heating region. 


.15) 

.16) 

.17) 
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2.7,2 Thermal Balance in Shaper Coil 

In this section* we consider the problem of heat production and 
dissipation in the fine shaping coil, which we will consider to be a simple 
hairpin. We will show that wires small enough to put close enough to the 
melt for effective shaping will not be able to lose enough heat by radiation 
to keep from melting at the current densities required. We will also con- 
sider the possibility of cooling the wire by making it a hollow tube with 
fluid flowing through it, and show that there is some possibility of cooling 
with flowing gas. We were not able to make a realistic estimate for cooling 
with a liquid flowing through the tube, but believe it might be feasible. 
Cooling might also be accomplished by partially embedding the coil in a 
high- thermal- conductivity , low-electrical-conductivity heat-pipe material; 
this possibility has not been investigated in any detail. 

We will make the following assumptions: 

1) We will neglect the large incident thermal flux from the ribbon itself, 
and only consider the heating by the current in the wire. Although in some 
circumstances a good part of the incident flux (peaked around 1 to 1.7 pm) 
might be reflected, in general the only justification for this assumption is 
that it makes the problem much easier to solve, and in a more complete 
analysis it ought to be included. 

2) The radial dependence of the ac current density can be ignored as far as 
the thermal, rather than electromagnetic, part of the analysis is considered. 
This can be shown to be a very small effect. 

3) The temperature dependence of the resistivity and thermal conduct iv: " 
of the shaper wire can be neglected. It car be seen that if these effec 
are not negligible, they will make the wire more difficult to cool. 

4) The wire radius a is much larger than the skin depth in the wire. The 
bend in the coil can be ignored. These are safe assumptions. 

5) The wire is made of copper, of nominal properties as follows: 
melting point, T - 1356 K, 

_ Q 

resistivity p^ c = 1.7 x 10 ohm-m, 
skin Hpnth ft a 14.3 um fat- 71.1 MHz), 

* CU ' 

thermal conductivity K- = 3,9 x 10^ W m” ^ K" ^ 

Cu 
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We will suppress the subscripts ("dc" and "Cu") in the equations. It is 
doubtful that a better metal than copper can be found for the wire. As 
frequently mentioned, a practical shaper may also include some sort of 
graphite or similar concentrator, whose thermal properties should be in- 
cluded in the overall analysis; but it is instructive to see just what is 
required of the coil itself. 

The steady- state temperature in the wire will be given by an 


equation of the form 


(25) 


T * A + B In r • 


A 0 r 


4K 


where r is the radial distance from the center of the wire and Aq is the 
power generated per unit volume. For a » 6, 

d 1 ”) 2 . . 


2 3 

2ttV 6 


K ill 

K dr ' r = a 


€ E (T 


C 


<3 


this equation is valid whether the wire is hollowed into a tube or not. 

A and B are constants to be determined from the boundary conditions. We 
will examine in turn the two limiting (and relatively simple) cases of 
cooling by radiation only and cooling by forced convection in a tube. 

I. Solid Wire Cooled by Radiation at Surface . In this case, B 
must be zero and A is determined from 

,4 

r * a^ * * 

where e is the emissivity of the copper surface, £ is the Stefan- Boltzmann 
constant, and we have assumed radiation into a medium at zero temperature. 

We readily obtain 

/ A o a \ 1/4 

+ VST » • < 2 


A , 2 2. 
A Q (a - r ) 


4K 


The r->ximum temperature at the center of the wire may be written as 


T = T. + T 
max 1 s 


(2 


whe re 


T l = A Q a /4K 


and T 


s 


the surface temperature , = (A 0 a/2e£) 


1/4 


.7.22) 

.7.23) 


2.7.24) 

.7.25) 

.7.26a) 


(2.7.26b) 

(2.7.26c) 
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It is not hard to show that for the situations of interest will be less 
than one degree and can safely be neglected. (To be explicit, 

T x - 3.83 x 10' 5 (i rms ) 2 /a 

where is in K, i mS in A , and a in ram) . Also we find that for a shiny 
copper surface (e «^.02), the wire will immediately fuse at any power levels 
we have considered; so we will assume the surface is oxidized (e 6) or 
covered with some highly emitting material (e ^ 1), This of course will 
make the surface a better absorber of radiation from the molten Si, but as 
previously stated, we will not attempt to take this into account. By setting 
T g * T , we can calculate roughly the minimum wire radius a for a given 
current. The current, in turn, can be related to the wire-to-melt distance 
b for shaping the ribbon to 400 um thickness at 21.1 liiz, the conditions 
principally discussed in earlier sections. The results are shown in 
Fig. 2.7.2. Of course unless a < b, the wire cannot be fitted in next to 
the melt. It is apparent that even for e « 1, the wire must be made quite 
large and moved rather far away from the melt if it is not to melt itself. 

The overall power requirements for the various situations were discussed 
in Sec. 2.5. Since ail our assumptions have been such as to neglect effects 
that would tend to raise the temperature more than we have calculated, it 
appears that it will not be possible to rely on radiative cooling of the 
shaper wire. 

II. Hollow Tube with Internal Fluid Cooling . We assume the 
shaper is a cylindrical tube of interior radius a^ and exterior radius a. 

We will now assume there is negligible radiation at the outer surface (e % 0) , 
to simplify the discussion. At the inner surface, we will take as the 
boundary condition 

dT, h / T \ (2.7.27a) 

dr ' r - a x ‘ K V r = 3i f ; * 

where the heat transfer coefficient 

h = K^Nu/2a^ , (2.7.27b) 

being the thermal conductivity of the cooling fluid and Nu being the 
Nusselt number, a dimensionless parameter which we discuss below. T^ can 


i 
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FIG. 2.7 


b ( mm) 


.2. RADIATIVE COOLING OF SHAPING COIL. MINIMUM COIL RADIUS 
TO AVOID FUSION, a, PLOTTED AGAINST COIL-CENTER-TO-MELT- 
SURFACE DISTANCE, b, FOR SHAPING RIBBON 0.4 m THICK AT 
21.1 MHz. « - EMISSIVITY OF COIL SURFACE. 
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be taken as the average temperature of the cooling fluid, provided the 
outlet absolute temperature is not more than 3 times that at the inlet 
Evaluating A and B from the boundary conditions, we find 



+ In 1 



(2.7.28) 


From this equation, we can show that the temperature difference across 
the wall, 


T(r»a) - T(r»a^) < T^, where T^, defined by (2.7.26b), is less than 

one degree. Thus we consider only the temperature at the inner wall, and 
write 

T(r- V - T f • AT - Tl fS_ [l -£) 2 ] . P.).!)) 

If the temperature variation of the viscosity is not too strong, the 

/oft 07\ 

Nusselt number for most gasses is well approximated by*' * 

DO 

Nu = .013 Re , (2.7.30) 

3 

in the turbulent flow regime (Reynolds numbers Re > 3 x 10 ). The 
Reynolds number is given by 

Re = 2 a] v/u , (2.7.31) 

where v is the average flow rate of the gas and u is its viscosity at the 
average temperature. Defining the ratio X = a^/a, we see we can write 

AT _ 2’ 17 K / u V 83 1-X 2 . (2.7.32) 

T, “ .013 K f \av / x .83 

X must be chosen to make the shaper sufficiently strong; let us assume 
a ■ 1 mm and X * .5. Also let us assume that the gas is air, precooled so 

4 

that the average temperature is room temperature (294K) . Then K/K* fe 1.6 x 10 

2 r 
and u .16 cm /sec; so 

AT 3.84 x 10 4 
T 1 * V 83 

where v is given in m/sec. For turbulent flow, v 

3 

about 30 m/sec. At that rate, AT/T^ * 2.3 x 10 . 


will have to be at least 
For a * 1 ran and b a 1.4 mm, 
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as .46 K, lienee the wall temperature would rise to 1336 K, or just a bit 
below the melting point. Since this might be brought down by higher flow 
rates and by conduction and radiation losses, it seems as if cooling with 
a gas might be possible for a shaping wire of the dimensions assumed, it 
is not likely that smaller wires closer to the melt could be cooled this 
way though, despite their carrying smaller currents. We attempted to make 
a similar analysis for cooling by a liquid, but found ourselves in an 
Inconsistent range of parameters for use of the appropriate generalization 
of (2.7.30). One would expect, though, that if cooling with a gas is 
feasible, cooling with a liquid at lower flow rates would also work. 



2.8 Summary and Conclusion 


We have described methods of calculating fields and forces 
induced in a silicon ribbon by a simple hairpin-shaped coil. We have 
determined the current and power necessary to produce ribbon of a speci- 
fied thickness by electromagnetic shaping, both for an isolated ribbon 
and for a ribbon being drawn from a nearby levitated spherical melt In 
both cases the current required for shaping to .4 mm thickness is larger 
than has been previously thought-- in the range of 100 A. For the isolated 
ribbon, the power required is also large, but for the ribbon drawn from 
the levitated sphere, the total power is only a few kilowatts. Departure 
of the melt from a flat ribbon shape, location of the freezing interface 
far from the shaping plane, and excessive heating of the ribbon by the 
eddy currents may all be severe problems which will have to be overcome 
in any practical system. Numerous simplifications and approximations have 
been made in the course of these calculations, but it is extremely doubtful 
that more accurate calculations would alter qualitatively any of the con- 
clusions presented here. 


3. OTHER INNOVATIVE PROCEDURES 


3.1 Introduction 

The object of this task ..as been the evaluation of possible means 
for shaping silicon ribbon other than using electromagnetic forces, and making 
a detailed analysis of those methods which appear most feasible from the stand- 
point of space production. Emphasis was placed more on developing understand- 
ing of the physical mechanisms of importance and on establishing high- feasi- 
bility concepts for shaping than on engineering design. 

By way of introduction, we shall describe three techniques used for 
shaping materials (terrestrially) from the melt: the web-dendrite, the EFG 
(Edge-Defined Film-Fed Growth), and the Pi Iking ton Float Glass methods. 

Web Dendrites 


In this process, a seed crystal is induced to form two dendrite arms, 
which are curved down into the melt by surface forces imposed by the liquid 
material. Soon after the start of growth, the two dendrites form what is 
essentially a parallel frame, supporting a web of liquid which freezes shortly 
above the surface of the liquid. The web material has a twinned structure, 
and the crystallographic orientation in the plane of the sheet is off a close- 
packed orientation as a consequence of the dendritic growth process. Growth 
rates are of the order of 5 cm/sec. One of the problems associated with this 
method is that, as a result of liquid drainage induced by gravity, the material 
does not have a uniform thickness. Hence, the cooling rate is not uniform, 
and thermal stresses are induced unless cooling rates are extremely well con- 
trolled. Generally, at web widths of about 3 cm the dislocation structure 
induced by thermal stresses is so great that the web becomes polycrystalline. 
Furthermore, the dendritic growth process results in variations in the con- 
centration of alloying elements across the width of the web. In addition, 
bet jusc dendritic growth is uncons t rained , the process is rather delicate to 
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EFG 

Liquid material is drawn into a die by capillarity forces. From the 
pool in the top of the die, it is drawn out by a seed crystal in a manner 
analogous to the Czochralski method. The liquid departs from the die with 
a shape defined by bounding vertical (in the direction of pulling) surfaces. 

In principle, if this method were developed properly, the material 
would have properties closely comparable with the best Czochralski material. 

In this case, it would suffer from two main drawbacks. First, the amount 
of yield would be limited by stability considerations. This drawback is over- 
come, to some degree, by using multiple grcwth slots. Second, the dies wear 
excessively, and the die cost forms a significant amount of the material cost. 

We see another possible disadvantage for future development, l.e., there 
is a lower limit to the ribbon thickness produced by this method, which we 
estimate to be about 120 urn. 

Filkingtoa Float Glass Process 

Molten glass flows out of the melting tank, over the surface of a 
tank of molten tin. As it flows over the tin its temperature gradually falls 
until it becomes solid. It then passes over a weir, onto a conveyer, where 
It is cut into convenient lengths. 

3.2 Factors Affecting Ribbon Formation 

Let us now consider the relative importance of the various factors 
which affect ribbon formation in space. In particular, we Include (1) 
capillarity effects, which tend to alter the shape of the liquid until mini- 
mization of net interfacial energy of the entire system is achieved, (2) viscous 
effects which tend to retard the flow of the liq '.d, (3) fit** forces resulting 
from the fact that the liquid is in motion, and (4) accelerational effects with- 
in the liquid. We shall demonstrate that, for liquid silicon, the first of 
these factors is dominant over the last three. 


i 





3,2.1 Capillarity Forces Versus Viscous Forces 

Let us first compere the relative influence of capillarity forces 

in a liquid-silicon system, as opposed by viscous forces. We shall consider 

a relatively simple example for which a quantitative treatment Is available, 

l.e., the spreading of a viscous liquid drop on a planar solid surface under 

the action of Interfacial forces. We use here the results of the theory of 

(28) 

"vetting” kinetics developed by Strella v ' and extended by Markvorth and 

(29) 

G leaser . They considered the problem of a liquid drop, initially 
spherical and touching the planar surface at only one point, as it subsequetly 
spreads over the surface, finally attaining its equilibrium configuration. In 
particular, they shoved that 


4 1/3 T) R 


f(e,ej 


where Of is the surface tension of the drop, T | its viscosity, R q its initial 
radius, t is the time measured relative to the initial condition, for which 
t » 0. 6 is the Instantaneous contact angle between the drop and the surface 
and is measured relative to the planar surface not in contact with the drop, 
and 0^ is ^he equilibrium contact angle. (The contact angle 0 approaches 
0^ in the limit as t approaches infinity it can be seen, from Eq. (3.2.1), 
that, as expected, high values of or lead to correspondingly high rates of 

approach to equilibrium, whereas high T| leads to slow rates. The functional 

(29) 

form of 1(0,0^) is unimportant for the present work but one could shew, 
with the help of results reported in Refs. (28) and (29), that the spreading 
process is generally largely complete at times such that this function has 
a value of roughly 5. Hence, defining < as the time required for the function 
to attain this value, we find, 


4 l/3 T| R /a 
o 
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Using the values R q ^ 1 cm, 7] ar 0.02 po*:.®, and a 720 dyne/cm, we find that 
r Is of the order of 10~^ second, so that surface forces in ’ <quid silicon 
generally overcome retarding viscous forces extremely quickly. Thai this Is 
so is, however, not surprising in view of the relatively high su^ace tension 
of liquid uilicon and its relatively low viscosity — for example, the viscosity 
of water is only about a factor of two smaller than the value used here for 
liquid silicon. 

3.2,2 Capillarity Forces Versus Accelerations! and Flow Force s 

We shall use a somewhat different approach to evaluate tho* com- 
parative effects of capillarity forces versus accelerational and flow forces. 

In particular, we consider two of the various dimensionless parameters that 
are commonly used to identify behavior similarities amoug system* of different 
size and also to define different regimes of behavior; these tv o h ^rameters 
are the Bond number, Bo, and the Weber number, We, defined as 

Bo s paL 2 /-y 

We £ pL V 2 /y 

where p is the fluid density, a its acceleration, V Its velocity, L a 
characteristic length, and y the fluid surface tension. The Bond number 
defines regimes in which either acceleration or capillarity forces dominate, 
whereas the Weber number compares flow forces and capillarity forces. The 

appropriate hydrodynamic regimes can be illustrated graphically (see Reynold.-*^ 
(31) 

and Otto v on a plot of Bo versus We, as shown in Fig. 3.2.1. Assuming, 

for liquid silicon, order-of-magnitude -estimated values for o, a, L, V, and 
3 2-3 

y of 1 gm/cm , 1 m/sec 10 g N , 1 cm, 10 cm/sec (roughly the anticipated 

3 -7 -3 

pull rate) , and 10 dyne/cm, rtspectively , we find We ~ 10 and Bo - 10 

The point in Fig. 3.2.1 characterized by these numbers can be seen to lie deeply 

within the “cap^ 1 lary-dominated" region, indicating again that capillarity 

forces in low gravity play a dominant role in controlling the instantaneous 

shape of the liquid-silicon system. 



77 



FIG. 3.2.1. HYDRODYNAMIC REGIMES (AFTER REYNOLDS AND OTTO). 1. lE QUANTITY 
Fr IS THE FROUDE NUMBER, WHICH IS EQUAL *0 (We/Bo ) . 


ORIGINAL PA*r 31 
OF POOR QUAIira 
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3.3 Use of Capillarity Forces to 
Shape Silicon Ribbon 


3.3.1 Capillarity-Induced Interaction Between 
Solids and Liquids 

From a very general point of view, the shaping of liquid silicon 

into ribbon, via capillarity forces, would be achieved through the interaction 

of the liquid with an appropriate arrangement of solid surfaces. Assuming 

that inter facial forces control the shape of the liquid, the shape would be 

that for which the net interfacial energy of the system (liquid plus solids 

with which the liquid comes into contact) is minimized. Unfortunately, the 

problems that one encounters in determining the minimum- energy configuration 

for a relatively complex system are generally exceedingly difficult. Indeed, 

(32) 

it has beer pointed out by Nitsche that only in 1951 was it proven that 
a sphere- like surface (certainly a simple system from the point of view of 
the more complex arrangements needed for ribbon shaping) having constant 
mean curvature must in fact be a sphere. To illustrate the approach that 
one can take in calculating capillarity-driven equilibrium solid-liquid 
configurations, we shaM consider here the simple, but important and pertinent, 
problem of calculating the equilibrium configuration of a quantity of liquid 
situated on a rigid, spherical solid surface. 

The system we are considering is illustrated schematically in Fig. 
3.3.1. Thi droplet shape, along the liquid-vapor surface, is characterized 
by constant curvature, and the locus of points describing the boundary between 
the solid-vapor surface and the solid-liquid interface is a circle. All 
angles and radii of curvature are taken as positive for the specific con- 
figuration il ,trated in Fig. 3.3.1, and their signs would be altered in 
appropriate fashion for concave interfaces. Let u3 suppose that the radius 
of the liquid system in a free, and therefore spherical, configuration is 
r. Then, it is relatively easy to show that the change of interfacial energy, 
AE^, brought about by "attachment 1 ’ of the liquid to the surface, is given by 
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AE t - 2np y^ v t 1 - cos (9 + or)] 

+ 2ttR 2 y\ ls - y sv j (1 - cos o) - 4nr 2 y ^ (3.3.1) 

where p is the radius of curvature of the liquid-vapor interface after 
attachment to the solid has occurred, R is the iadius of curvature of the solid, 
and y QV > and are the re. >ective free energies per unit area associated 

with the liquid-vapor, the solid-vapor, and the solid- liquid interface? . The 
angles 0 and a are defined in Fig. 3.3.1, from which we note that 6 is the 
’’contact angle" between the liquid and the solid. The angles are all taken 
as positive r or the configuration illustrated in the figure. 

Now the equilibrium configuration of the system is that for which 
AEj is minimu l (we are assuming zero gravity) subject to the constraint that 
the volume of the liquid remains constant upon attachment to the surface. 

To evaluate the nature of this configuration, we employ Lagrange's method 
of undetermined multipliers and define a quantity Y such that 

Y = AE. + V V ^ VJ (3.3.2) 

1 v a 


where X is a Lagrange multiplier, is the volume of the attached liquid, 
and is the volume of the free liquid. Now, one can readily show that 


V 

a 



e (e + o) 



(3.3.3) 


V 


f 



(3.3.4) 


where 


F(x) = 1/2 (1 - | cos 


+ 1/2 cos x) . 


i 
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la addition, we see from Fig. 3.3.1 that 


R sin a ■ p sin (8 + a) 


(3.3.5) 


Combining Eqs.(3.3.3) and (3.3.5), 
4nR 3 • sin 3 ar 


1 = * ! ? (6 + a) - ^(a) 

a 3 > sln J (e+a) * 


Next, combining Eqs.(3.3.1) and (3.3.5), 

2 Y iv 8ln2 * 

AE - 2ttR [1 - cos (0 + a) ] 

sin (9+0) 


(3.3.6) 


Kyiv - w (1 - C08 “>j - 4nr Y 


tv 


(3.3.7) 


Eqs . (3.3.2), (3.3.4), (3.3.6), and (3.3.7) can be combined to yield an 
expression for Y entirely in terms of the variables ot and 0 (treating r, R, 
and the interfacial energies as fixed parameters) . 

The minimum value of AE^ is determined by setting 

d * = I de + S da = 0 * (3 - 3 - 8) 

and hence by equating each of the partial derivatives in (3.3.8) to zero. 

First, setting dY/<30 ■» 0 and solving for X, we obtain 


X - sin (0 + ap 

R sin ot 


(3.3.9) 


It is interestirr here to observe that Eos. (3.3.3) and (3.3.9) can be com- 
b’ned to yield 

X = - 2y'p 

so that X is simply the negative of the "Laplace pressure” associated with 
the attached liquid. Next, we set ?)Y/*3r = 0 and subsequently eliminate 
X in the resulting expression by means of (3.3.9). We thus obtain 


v = V, + v. cos 8 
Y sv T /s tv 


(3.3.10) 
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Eq. (3.3.10) which is identical to the familiar equation of Young and 

Dupre ', describes the relation between the equilibrium contact angle and the 

various interfacial energies. This expression has been derived here for the 

special case in which the solid surface is spherical, although the same form 

for this relationship has been obtained for other situations as well (e.g., 

( 33) 

liquids contained within closed axisyimnetric containers v ') . It will prove 
to be essential for use in evaluating the equilibrium configuration associated 
with the more complex systems which would be used to shape liquid silicon 
into ribbon. 

We note, in concluding this analysis, that Darbro^’^ has 
discussed a number of important concepts relative to phenomena associated 
with the nature (particularly under low gravity conditions) of such things 
as liquids, liquid films, liquid-solid interactions, and the solidification 
of liquids. 

3.3.2 Melt Containment 

For the shaping of ribbon from a bulk melt, it will be important 
to control the shape of the liquid- vapor interface at the melt surface. 

Because of the dominance of capillarity forces, the shape of the interface 
is controlled by the shape of the container as well a a the contact angle 0. 

To illustrate the manner in which shape control can be achieved, we consider 
the problem of obtaining a flat liquid-vapor interface, although the procedures 
we describe apply equally well to obtaining interfaces having non-zero curva- 
ture. We shall assume that the Young-Dupre' equation, as derived above for 

the interaction of a liquid with a spherical (concave or convex) solid surface, 

(33) 

applies to other curved surfaces as well 

One possible configuration, for which the container is spherical, 
is illustrated in Fig. 3.3.2. The principal drawback here Is that a flat 
liquid-vapor -aterface can be obtained for only one depth of liquid, as 
illustrated in the figure. If the container holds either more or less then 
this particular amount of liquid, then the liquid-vapor interface will be 
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curved (concave upward for the former case , convex upward for the latter) . 

A more practical arrangement is shown in Fig. 3.3.3. Here the container is 
a surface of revolution, the walls forming a circular cone over at least 
some region, the apex angle 4 of the cone, as illustrated in the figure, 
being n - 20. For this situation, the liquid-vapor interface will be flat 
for any liquid depth as long as the interface intersects the container with- 
in the region that is circularly conical. 

The above examples illustrate the importance of container shape 
in obtaining liquid-vapor interfaces having given desired curvature. It is 
clear, however, that the problem is one which can be met by suitable choice 
of container material and shape. 


3.3.3 Capillarity-Driven Ribbon Formation 

We shall now describe a general concept for ribbon formation based 
on the assumption that interfacial-energy effects control the shaping process. 

In particular, we shall utilize a crystallization process which is "constrained"*, 
so that the material produced can, in principle, be of the highest quality. 

The two key features of the process are the following: 

1. A framework is used to draw liquid sheet from a bath in a manner 
analogous to the role played by the dendrites in the web process. 

2. Constrained growth is coupled with rapid rate of production by pulling 
the crystal in a direction aligned at an angle to the growth direction. 

(We note that the benefits offered by the low gravity environment of an 
orbiting spacecraft relative to constrained crystal growth have been discussed 
elsewhere^ for the particular case of the growth of beta alumina from the 
melt under low-g conditions.) 


* Constrained crystal growth is tha* , >r which the molten material ahead 
of the solid-melt interface is *■ . e, i.e., above the melting p^int, as 

opposed to the case of unconstra * * crystal growth for which the liquid 

ahead of the growing crystal is metastable, i.e., below its equilibrium 
freezing temperature. Freezing occurs in the constrained- growth situation 
as a result of the isotherm corresponding to the freezing temperature being 
moved relative to the melt. Maintenance of constrained growth prevents the 
solid-melt interface rrora breaking down. If, however, the growth rate be- 
comes too fast, the interface becomes cellular; and in the limit of even 
faster growth rates, it becomes dendritic. If, however, the growth rate 
becomes too fast, the interface becomes cellular; and in the limit of even 
faster growth rates, it becomes dendritic. 
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FIG. 3. 3.3. RELATION BETWEEN APEX ANGLE 0 AND CONTACT ANGLE 6 IN ORDER 
TO ACHIEVE A FLAT LIQUID- VAPOR INTERFACE INSIDE A CONTAINER 
WHICH IS A RIGHT CIRCULAR CONE OVER A PORTION OF ITS SURFACE. 
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Aspects of this process are discussed belov in more detail together 
with our view of some of the attendant problems. The concept of "off-angle” 
growth, mentioned under (2), above, is sufficiently important in its own 
right to be described in a separate section, which follows below. 

The basic procedure for utilizing capillarity forces to produce 
ribbon is illustrated schematically in Fig. 3.3.4. The melt is drawn up in the 
desired ribbon shape, along a frame which defines its boundaries. The frame 
could, in principle, either be stationary or move with the liquid at a rate 
such that the liquid ribbon does not move relative to the frame. A stationary 
frame would result In a simpler experimental apparatus, but a movable frame 
would undergo less wear due to frictional effects resulting from the surrounding 
moving liquid; in addition, deleterious bound ary- layer effects in the moving 
liquid would be minimized if the frame is also moving. Solidification of 
the liquid could be induced either before or after the ribbon is separated 
from the frame. Some particular aspects of this concept that would require 
further consideration are described below: 

Film Stability . Viscous forces, as we have shown, tend to be 
small compared to capillarity forces; but they do exist and do tend to oppose 
film formation, and the extent to which these forces constitute a threat to 
film stability constitutes a problem which must be analyzed. For the web 
dendrite process, which in this respect is a close analog of the proposed 
process, pull rates of about 5 cm/sec have been attained despite the super- 
cooled (and thus more viscous) melt. We anticipate, for the space process, 
that pull rates of the order of ,01 to .02 cm/sec will be attained. 

A second stability question is involved with the stability of a 
thin film or sheet. Without externally applied forces, the liquid film 
supported by a framework should be stable. If subjected to external 
vibrations, it could support limited amounts of longitudinal and/or transverse 
standing-wave systems; these waves would be damped, of course, by viscous 
flow of material. We did not consider this problem, but the limits of applied 
vibrations, beyond which the liquid sheet becomes unstable, muse eventually 
be determined as a function of sheet dimensions. The success of the web 
dendrite method gives us some confidence that the potential stability problems 
associated with our process can be satisfactorily resolved. 
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Materials for Frame and Crucible , Materials developed for the EFG 
process should be directly applicable to our work. The most successful has 
been graphite, which seems to be an obvious first-choice material. Other 
materials which have been successful are silica and silicon carbide. These 
materials may need surface treatment in order that they may be "wetted" by 
the liquid silicon. 

Mechanical Stresses . The greatest concern regarding the concept 
ve have developed is related to the mismatch in thermal -expans ion properties 
between the frame and the single crystal formed within it. Our concerns on 
this ground have been reinforced by the experience of the web-dendrite crystal 
growers. Our response to this problem is twofold: First, an examination is 

called for of the difference in the stress distribution be^vaen a very thin, 
wide, film and the less flexible thick, narrow web crystal. Second, some 
potential practical solutions to this problem exist. These solutions include 
cutting away the frame closely above the solid-liquid interface region by 
some mechanical procedure, and the alternative in which the frame is drawn 
away Just before the freezing interface region. The latter method would 
provide a less flat sheet than one formed between a frame, but would com- 
pletely remove the thermal- stress problem. 

Provided that the frame-stress problems can be sacisfac torily 
removed, other potential sources of stress jeem to be less important. For 
example, since the web crystal growers coiled their crystals continuously or 
a 90 cm reel, the crystals formed in the space process, which should be 
thinner, may be coiled at least as tightly. 
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3 .4 Off- Angle Growth 


The concept of off- angle growth, mentioned briefly above, offers 
some important advantages for the production of single-crystal ribbon in 
space. Suppose a single-crystal ribbon is oelng grown such that the pull 
velocity and growth velocity are aligned in different directions. Let V 

P 

be the pull velocity, V the growth velocity, and 0 the angle between the 

© 

pull and growth directions. (The quantity 0 defined here is distinct from 
the "contact angle" discussed above, which was also denoted by 6). This 
concept is illustrated schematically in Fig. 3,4.1. For the configuration 
shown here, as well as for those illustrated in Figs. 3.4.2 and 3.4.4, it is 
assumed that the pull direction lies in the plane of the molten sheet; this 
need not be the case, however, as we shall discuss below (and as is shown 
in Fig. 3.4.3). 

We can use Fig. 3.4.2 to help derive the relationship between 

V . V , ar.d 0. Let us define A as the actual area of the solid-melt inter- 
P' g g 

face (i.e., tb*t normal to the growth direction) and A^ as the interface 
area projected normal to V . In order that mass be conserved, we require 
that 

V A = V A 
8 g P P 

But , one can readily see that 


so that 


A = A cos G 
P g 


V = V cos 0 
g P 


Hence, non- zero * .iue of 9 would result in an enhancement of pull rate, 


i.e., V 


V . 
g 


Two possible approaches to achieving off- angle growth are snown 
in Fig. 3.'. In one of these, the growth velocity vector is aligned at 
an angle to the plane of the molten ribbon; in the other, it lies within 
this plane. (Fig. 3.4.1 and 3.^.2 correspond tt the former case.) The 
latter situation is probably simpler to achieve experimentally, since heat 
flow within the molten ribbon would occur along the plane of the ribbon. 


I 






FIG. 3. A. 2. ILLUSTRATION OF QUANTITIES USED TO DETERMINE RELATIONSHIP 
BETWEEN V„. Vp, AND 6. AN EDGE VIEW OF THE MOLTEN RIBBON 
IS SHOWN TOGETHER WITH THE SOLID-MELT INTERFACE. 
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Other problems would exist, however, such as those resulting from the 
manner in which heat is radiated from various portions of the system; 
these would have to be treated in subsequent studies. 

Advantages to be gained from utilizing off- angle growth for 
ribbon production include the following: 

(1) Enhancement of pull rate. 

(2) More homogeneous crystal. Defects (e.g., spiral dislocations) tend 
to propagate along the growth direction and hence would terminate at the 
ribbon surface for off-angle growth. 

(3) Thinner ribbon. For the case illustrated in Fig. 3.4.2, for which 
the solidified ribbon and molten ribbon are not coplanar, the solidified 
material is thinner than the molten ribbon. 

(4) Well suited to space processing. The off-angle growth process would 
be difficult to maintain stably on earth. 

The portion of the solidified ribbon which, at the melt- solid 
interface, is most pointed, has some of the properties of a growing 
dendrite, and existing dendritic- growth theory is applicable. A sharply 
pointed surface would be unstable, because of the Gibbs-Thomson effect, 
and it would tend to become rounded as illustrated in an exaggerated 
manner in Fig. 3.4.4. The degree of rounding depends, in some measure, 
upon such factors as the rate of advance of the crystal and upon the 
crystal-melt inter facial energy. We expect the radius of curvature to lie 
in the range 1 to 5 pim. One would normally expect an interface of this 
nature to become irregular; however, thermal constraints would prevent this 
from happening in the present situation. 

It should also be noted that our concept of off- angle growth 
would lead to a surface exposed along the plane of the sheet which may not 
be a close-packed orientation. Such an orientation is obtained from web 
material . 

Off- angle growth is thus a potentially useful means for enhancing 
the rate of ribbon production in space as well as the homogeneity of the 
material. It should be emphasized, however, that this process is not 
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essential to the concept described above involving the use of capillarity 
forces for ribbon shaping. It would be used only if found to be feasible 
within available technology. 


3.5 Conclusion 


The processes described above, i.e., the use of capillarity forces 
to shape silicon ribbon and the application of off-angle growth to produce 
ribbon at faster rates and greater degree of homogeneity, appear to be 
feasible from a conceptual point of view. Engineering design problems do 
exist, such as maintenance of a stable configuration during off-angle 
growth, and these must be considered as part of further developments. 
Certainly, the low- gravity environment would have considerable influence 
upon design considerations. 

Additional shaping concepts that also appear to be feasible need 
to be evaluated against the processes proposed herein. (Indeed, we have 
considered a number of possible concepts, subsequently discarding them 
in favor of those described above.) For example, one potential process 
would involve zone melting, utilizing a spatially fixed molten zone. Poly- 
crystalline sheet would be fed along a given direction, melted within the 
zone, and single crystal ribbon pulled out from the molten zone. This 
concept could lead to problems, however, in obtaining a flat solidified 
sheet. It should be noted that off-angle growth could be utilized with 
this type of procedure. 
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4. CONCLUSIONS AND RECOMMENDATIONS 


As part of our overall assessment of the prospects for space 
production of silicon single crystal ribbon, ve reviewed the recent liter- 
ature relating to potential large-scale use of silicon for solar cells and 
had discussions with a number of researchers active in this field. The 
results of this survey have been considered In arriving at the following 
conclusions : 

1. There is no doubt that the need exists for much better methods of 
producing silicon ribbon for solar cells, if photovoltaic conversion is 
ever to emerge as a viable means of producing large amounts of electricity. 
Not only are much lower costs required, at almost every stage of the process, 
but also minimally low energy expenditures, in comparison to the total energy 
produced by the material over its useful life. 

2. There exist a number of competing technologies for production of Si 
ribbon. Some of these were described in Section 3.1, and others are under 
development. No process or concept at present has demonstrated any clear 
super lor i ty , though . 

3. No material other than silicon is given serious consideration for large- 
scale terrestrial solar power conversion, because of 1) availability of raw 
material, ii) good conversion efficiency, and ill) overall knowledge of and 
experience with the material. 

4. As far as production of silicon wafers for integrated circuits goes, it 
appears that any new process will have to be Justified strictly on economic 
grounds, and moreover will have to dislodge an entrenched, if rather awkward, 
technology. 

5. Current requirements for electromagnetic shaping of silicon ribbon are 
substantial even if small shaping colls very close to the melt are used 
Currents in the 100 A range will be required for shaping to 0.4 inn thickness 

6. Power requirements for drawing an electromagnetic ally shaped ribbon from 
a levitated spherical melt are moderate— a few kW, Including losses in both 
the shaping coll and in the mo] ten silicon, if the shaping coil is small. 
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7. Power requirements for thinning to a few hundred microns a ribbon which 
has been previously drawn by some other process are high. Such processes 
do not at present appear feasible with the amount of shaping power that is 
likely to be available in a space processing facility. 

8. It cannot be expected that electromagnetically shaped ribbon will be of 
uniform thickness. The thickness near the ribbon center could easily be twice 
that near the edges. This might not lead to major problems in applications, 
but could be a problem in keeping the ribbon from touching the shaping coils. 

9. Shaping coils could be cooled satisfactorily by convective transfer to 
fluid flowing through a hollowed wire, and possibly by other means. 

10. Eddy-current heating of the molten silicon will raise Its temperature 
significantly and will move the liquid-solid Interface rather far from the 
region of maximum shaping force, under typical shaping conditions, unless 
some form of substantial auxiliary cooling of the ribbon can be provided 

11. In the absence of electromagnetic forces, the principal mechanism 
determining the shape of liquid silicon in the space environment is capillarity, 
which dominates over viscous, accelerational, or flaw forces. The morphology 

of the liquid silicon is thus principally dependent on the shape of any 
solids with which it comes into contact and the equilibrium contact angles 
between the silicon and the materials of these solids. Ribbon could in 
principle be produced simply by drawing between parallel wires. 

12. The rate of ribbon production, the homogeneity of the solid material 
produced, and the thinness of the shaped ribbon can all be enhanced by a 
process of off-angle growth in which the ribbon Is pulled at an angle to the 
direction of solidification. 

On the basis of these findings, we make the following recommendations : 
l. Although a number of possible difficulties in and drawbacks to space 
production of silicon single-crystal ribbon have been identified, nothing 
has been found which definitely nullifies the technical feasibility of such a 
process. Therefore we do not reconxnend that work in this area be halted. 

In view of the overall importance of being able to produce high-quality silicon 
ribbon at reasonable expense, i recomnend that work along these lines con- 
tinue to the extent justifieo by the overall goals of the space-processing 
program. 
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2. The concepts of greatest promise, both jith regard to electromagnetic 
shaping and to other methods, are very "space-oriented" and there is no 
reason to believe they can be convincingly demonstrated or modeled on earth. 
Obviously, therefore, terrestrial experiments cannot be generally recommended, 
although some experiments to test a particular point might be desirable. For 
example, it appears that a realistic investigation and test of many aspects 

of the capillarity-driven shaping concept described in Sec. 3.3 could be 

carried out on earth using immiscible liquids of similar densities, as in 

(32) 

the classic experiments of Plateau. 

3. A number of specific problems associated both with electromagnetic and 
non-elec tromagne tic shaping of ribbon have been pointed out in the conclusions 
above and elsewhere in the report. The primary need at the present time is 
for Innovative conceptualization to develop practicable ways these problems 
might be dealt with. If and when such ideas are formulated, further modeling 
to obtain a more quantitative assessment of their validity would be called 
for. Only after development of a demonstrably sound procedure for ribbon 
shaping would further design of the overall space processing system appear 
appropriate. 
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